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Executive Summary 


Investigation and Feasibility Assessment of TOPAZ-II Derivative 

for 

Space Power Applications 
Using a Gain Scheduled Controller Design 

The use of shunt regulators for load-following of proposed static space nuclear power 
systems (SNPSs) such as TOPAZ raises a number of concerns, such as the possibility of a 
failure in the shunt regulators requiring reactor shutdown, or the possible need to deliver 
somewhat higher power level to the load than originally expected. Therefore, a back-up 
system is needed in TOPAZ to eliminate the possibility of a single-point failure in the shunt 
regulators, and to increase the overall system reliability despite changing mission needs and 
deteriorating equipment. 

The objective of this research is to develop a controller for the voltage regulation in static 
SNPSs, which is capable of overcoming system variations resulting from operation at 
different power levels. This is accomplished by developing several linearized models of 
TOPAZ using a simplified, validated reference model of an integrated nonlinear SNPS 
model, covering its entire operating envelope. A dynamic compensator is designed using 
each of the linearized models based on the Linear Quadratic Gaussian with Loop Transfer 
Recovery (LQG/LTR) method. The LQG/LTR method is a systematic multivariable control 
system technique with several desirable properties such as a nominal stability and 
performance, as well as stability robustness to unmodeled dynamics. The various 
compensators matrices, including the gain and system matrices, are fitted to a scheduling 
variable, namely the SNPS electric power produced, to obtain a nonlinear gain-scheduled 
compensator. 

The performance of the gain-scheduled compensator is systematically investigated via 
transient and steady-state simulations using the integrated nonlinear SNPS model. The 



simulations, which include sensor noise at the plant output, demonstrate the effects of 
variations in the fuel temperature reactivity feedback coefficient on the load-following 
capabilities of the SNPS. Robustness analysis results of the gain-scheduled compensator 
demonstrate that the proposed control concept exhibits a significant degree of operational 
flexibility, and it is primarily intended for long-term space mission requiring significant 


levels of power system autonomy. 
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CHAPTER I 
INTRODUCTION 

The objective of this research is to develop a gain-scheduled adaptive control system 
that achieves voltage stability in a proposed static energy conversion space nuclear 
power system (SNPS), while generating different levels of electric power. The SNPS 
being analyzed in this research uses a fast, lithium-cooled nuclear reactor with thermal 
power output of 6.8 MW t and which has been proposed for use in future space 
mission by various researchers [1], [3]. This SNPS, capable of providing 300 k\V t 
electric power output, is based on the General Electric (GE) Company thermoelectric 
power conversion design for the SP-100 system [1]. 

The systematic controller design is based on linearized models of the SNPS valid 
at different equilibrium points. The linear models are generated from a more detailed 
nonlinear SNPS model, which has been previously utilized in the design of a model 
reference adaptive control algorithm. This controller enabled an integrated SNPS 
model to follow the predictable and desired response of the nonlinear reference model 
used in designing it, despite changes in the SNPS operating parameters throughout 
its entire lifetime [21. 

Linear controllers are designed using the relatively new method of control system 
design for Multiple Input- Multiple Output (MIMO) systems known as the Linear 
Quadratic Gaussian with Loop Transfer Recovery (LQG/LTR). The linear controllers 


This report follows a style based on the IEEE Transactions on Automatic Control. 
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are then scheduled (fitted) to arrive at a nonlinear controller with adaptive gains. The 

nonlinear controller is tested using the integrated nonlinear SNPS model throughout 
its operating regime. 

In this chapter, the motivation for the voltage controller design is discussed. 
A brief survey of previous work done on output voltage regulation of the proposed 
SNPSs is followed by the research objectives and contributions of this thesis. Finally, 
the organization of the remainder of this thesis is outlined. 

1.1 Motivation 

SNPSs have been under investigation as an alternative for both commercial and 
military space missions since the 1950’s [1], Such systems have a projected longevity 
of 7 to 10 years of continuous operation without the need for maintainance, they can 
supply a wide range of power (tens of k\V e to thousands of MW e ), and in certain 
power ranges they have significant mass advantage over other space power sources. 

In static SNPSs, the conventional approach for regulating the electric load 
demand has been the use of banks of shunt resistors (regulators). However, designs 
based on this approach raise a number of concerns, such as the possibility of a failure 
m the shunt regulators requiring reactor shutdown, or the possible need to deliver 
somewhat higher power level to the load than originally expected, because of changing 
mission needs or deteriorating equipment. Furthermore, operating the nuclear reactor 
at its nominal power level at all times, regardless of the load output of the system, 
results in the rejection of the excess load demand into space in the form of thermal 
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energy. This heat rejection increases the cooling requirements and consequently the 
total mass of the spacecraft. Results of parametric studies investigating the effects of 
multiple failures of the thermoelectric-electromagnetic (TE-EM) pumps and partial 
degradation of the waste heat rejection capability on the electric power output show 
that the system's power output is affected only slightly [13]. However, a failure in the 
shunt regulators has been pointed out to be a single-point failure. Hence, the use of 
shunt resistors in thermoelectric SNPSs calls for a back-up system. Such a back-up 
system could certainly improve the overall system reliability and eliminate a possible 
single-point failure. 

1.2 Literature Review 

Principal space reactor programs were terminated in the early 1970s, and only 
very limited activity took place during the remainder of the decade [1]. Around 
1979, the Department of Energy (DOE) initiated a modestly funded development 
program, focused on the heat-pipe-cooled reactor concept at the Los Alamos National 
Laboratory. By 1982, in view of the advent of a number of potential missions 
that would need nuclear systems, the DOE, the National Aeronautics and Space 
Administration (NASA), and the Department of Defense (DOD) had embarked on 
the development of a more active space reactor program [1;. Accordingly, in February 
of 1983, the SP-100 Program was initiated by NASA, DOE and the Defence Advanced 
Research Projects Agency (DARPA). The initial focus was on design of a 100 kW e 
system scalable from tens to hundreds of kW e . 
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Energy conversion in static SNPSs has been accomplished through thermoelec- 
tric or thermionic converters [6] , [7], [8j, [9], '10] . In SNPSs which employ thermionic 
energy conversion, control issues are similar to those using thermoelectric energy con- 
version [8] . Therefore a number of findings from this study are equally applicable to 
thermionic SNPSs. Thermionic energy conversion concepts have been applied to vari- 
ous space nuclear power systems under investigation, as alternatives to thermoelectric 
concepts [7], [8], [9], [10], 

The U.S. Air Force (USAF) has had a long standing interest in thermionic space 
power dating back to the early 1960s, when a heat-pipe-cooled thermionic converter 
was demonstrated through work at the predecessor to the Wright Laboratory (WL) 
[7], The STAR-C thermionic SNPS is a 40 kW e reactor system designed in response 
to evolving USAF space power requirements, which has eliminated the possibility 
of a single point failure resulting from a single converter failure [8]. The Idaho 
National Engineering Laboratory (INEL) has developed the Small Ex-core Heat 
Pipe Thermionic Reactor (SEHPTR) concep* to meet the safety and reliability 
requirements identified by the USAF [9], The SEHPTR concept has also eliminated 
the possibility of a single point failure resulting from a single converter failure. These 
studies indicate that converter failures can be tolerated in thermionic systems, because 
the converters can be interconnected in a series-parallel electrical network so that 
the current can flow around the failed cells through the remaining operational ones. 
Furthermore, the large number of independent power converters employed increases 
the network reliability with a small amount of additional redundancy [8], ;9j, 11]. 
Thus, high overall system reliability of the thermionic system has been achieved with 



0 


only modest converter reliability. Nevertheless, a potential single point failure due to 
shunt regulators is still present. 

The Advanced Thermionics Initiative (ATI) program was organized to inte- 
grate thermionic technology advances into a converter suitable for in-core reactor 
applications in the 10 to 40 kW e power range [71. The ATI program organization and 
its participating contractors have stated their goal for the development of a robust, 
long-life, high efficiency, simplified thermionic converter which can be incorporated 
into other space nuclear power reactor programs, and they have demonstrated their 
interaction during the annual program reviews. 

Since SNPSs, such as the thermoelectric-based SP-100. will frequently expe- 
rience changes in the external load demand, in 1985, a system simulation model. 
SNPSAM (Space Nuclear Power System Analysis Model), was developed at the Uni- 
versity of New Mexico’s Institute for Space Nuclear Power Studies, to investigate the 
load- following characteristics of an integrated SP-100 system without a shunt regu- 
lator for the electric load [4]. The results of this study showed that, although the 
nuclear reactor is always load-following, the integrated system is only load following 
up to a critical value of the external load demand; beyond this point, the system is 
nonload-following, requiring an external controller. 

In research investigating thermoelectric systems, the conventional approach for 
regulating the electric load demand has been to employ a shunt resistor (regulator) 
between the system and the payload [3], [4], |6], [12;, 13;, [14 . This approach, 
however, requires operating the nuclear reactor at its nominal power level at all 
times, regardless of the load output of the system, and the excess load demand is 
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rejected into space in the form of thermal energy [3'. This heat rejection, increases 
the cooling requirements and consequently the total mass of the spacecraft. In the 
reliability and vulnerability studies of the SP-100 system 13i, parametric analyses 
using SNPSAM-Mode-3 have been performed to investigate the effects of degradation 
of the system's waste heat rejection, due to a partial loss of surface area and/or the 
emissivity of the radiator panels, on the SP-100 operating parameters. The results of 
the parametric analyses assessing the effects of multiple failures of the TE-EM pumps 
and partial degradation of the waste heat rejection capability on the electric power 
output demonstrate that the system's electric output is affected only slightly 13 . 
However, a failure in the shunt regulators has been pointed out to be a single-point 
failure in these systems, and it has been concluded that the use of shunt regulators 
in thermoelectric, or even in thermionic, systems calls for a back-up load-following 
control system [3], [4], [13], 

A more recent study has made use of an adaptive controller to demonstrate that 
a model reference adaptive controller could cause a selected system state variable to 
track the transient trajectory of the corresponding state variable of a reference model 
with local stability in the SNPS [2], i5l. Although these results of the application 
of the adaptive controller suggest that a SNPS can operate in a load-following mode 
without the use of a shunt regulator, the particular application by Metzger et al. 2 . 
o does not address all of the questions that may arise from replacing a shunt regulator 
with a more automated, less hardware-intensive, load-following control system for a 
SNPS 3]. Such issues are, however, addressed in this study. 



1.3 Objectives of the Research 


The original objectives of this research are stated in the following sections. 

1.3.1 Linearized Model Development 

A number of linearized models will be generated by analytically linearizing a 
simplified nonlinear reference model of an integrated nonlinear SNPS model at various 
equilibrium points in an operating range of 50% and 150% of full power, along with 
numerically perturbing the latter model to obtain certain linearized reference model 
parameters at these equilibrium points. The linearized reference models generated 
are expected to be valid only for a finite range about the equilibrium points and 
will be validated with the integrated nonlinear SNPS model [17]. The number of 
linearized reference models to be generated will be determined by investigating the 
validity range of each linearized reference model and by assessing whether the entire 
operating range of the static SNPS has been covered. 

1.3.2 External Controller Design (Outer Loop) 

The external controller (compensator) designs for various operating power levels 
will be based on the corresponding linearized reference models of the plant and the 
LQG/LTR design methodology 16], [151. The LQG/LTR controllers are a subclass 
of the MBCs, which explicity incorporate the SNPS model dynamics into the design 
process. Hence, an accurate linearized reference model is required to design a high- 


performance external controller. 
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1.3-3 Gain-Scheduling of the Linear Compensators (Outer Loop) 

The linear controllers designed using the LQG/LTR design methodology are 
expected to be valid for operation within only a certain range about their respective 
equilibrium points. A single nonlinear controller (external controller) is therefore 
required, in order to continuously control the reactor in a range of 50% to 150% of 
the nominal power level. Gain-scheduling of the linear controllers, as well as of the 
linearized reference models, will be done by fitting each of the linear controller and 
linearized reference model parameters to a variable which is indicative of the system 
changing operating conditions. 

1.3.4 Design and Testing of the Adaptive Controller (Inner Loop) 

Model reference adaptive control is well developed for linear systems [2j, hence 
the linearized reference models of the SNPS will be used in the design of the adaptive 
controller. The constant gains in the reference adaptive control algorithm will be 
determined at various equilibrium points within the operating regime of the reactor. 
These gains will then be scheduled, to arrive at a single nonlinear adaptive controller 
that can cause the output (load voltage) of the SNPS to track the transient trajectory 
of the output of a reference model. The model reference adaptive controller will 
be tested with and without the external nonlinear controller, for the load-following 
of the integrated nonlinear SNPS model. As it turned out, the last objective was 
ne%'er implemented, because the outer loop of the dual-loop controller had satisfactory 
performance without the presence of an inner loop. 
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1.4 Contributions of the Research 

Following are the main contributions of this work : 

1) As opposed to current practice for load- following in proposed static SNPSs, a more 
automated, less hardware-intensive control system is proposed. This control system 
is intended as a back-up to shunt regulators, and is capable of effective control during 
variations in the power system response resulting from long-duration missions and 
component degradation. 

2) The designed linear controllers have certain guaranteed properties. Nonlinear 
transient response simulations indicate that to some extent these properties carry- 
over to the nonlinear gain-scheduled controller proposed in this study. 

3) The structure of an originally proposed dual-loop controller, comprising of an 
internal model reference adaptive controller and an external nonlinear controller, is 
proposed to be significantly modified. As a result of the performance indicators of 
the external nonlinear controller proposed in this study, there appears to be no need 
to still include the internal loop, i.e. the model reference adaptive controller, in the 
dual-loop controller for the type of parameter variations study in this thesis. 

15 Organization of the Thesis 

In Chapter II. the static SNPS under consideration is briefly described. The 
control problem in this particular system is stated. A dual-loop control mechanism 
is proposed as a solution for this problem, and current methods used to resolve the 
problem are presented. 
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In Chapter III, the various SNPS models used in this thesis are described. The 
generation of the linearized reference models from an integrated, nonlinear SNPS 
model is also described. The scheduling of the linear models is explained and 
evaluation of the scheduled model with respect to the integrated, nonlinear SNPS 
model is presented via simulation at various operating power levels. 

In Chapter IV. the linear control system design method used in this research is 
described and the mathematical preliminaries regarding this method are presented. 

In Chapter V. the design and gain-scheduling of the linear controllers proposed 
in this thesis is explained, and satisfaction of the design specifications by the linear 
controllers is verified. 

In Chapter VI, the performance of the linear and the gain-scheduled external 
controllers is analyzed at various operating power levels. The performance limits of 
the various controllers are determined, and a robustness study of the gain- scheduled 
controller is performed. 

In Chapter VII, the summary, the conclusions and some topics for further 


research are discussed. 



CHAPTER II 


VOLTAGE CONTROL IN STATIC 
SPACE NUCLEAR POWER SYSTEMS 


II. 1 Introduction 

In this chapter, the voltage regulation problem in static SNPSs is stated. A brief 
description of the static SNPS being analyzed is presented. The need for a back-up 
system for eliminating the possibility of a single-point failure in the shunt regulators 
of static SNPSs using thermoelectric, or even thermionic, power conversion is pointed 
out. A dual-loop control mechanism originally proposed to correct this problem is 
described. Furthermore, the conventional approach for regulating the electric load 
demand in static SNPSs and the studies of various researchers are briefly described. 

II. 2 Description of the Static Space Nuclear Power System 

The static SNPS under consideration uses a uranium nitride (UN)-fueled and 
lithium-cooled fast nuclear reactor with thermal power output of 6.8 This 

reactor system has been proposed for use in future space mission by previous re- 
searchers [1:, [3] . The particular system considered is capable of generating 300 k\V e 
and is based on the General Electric (GE) Company thermoelectric power conversion 
design for the SP-100 system [1], In this SNPS. energy conversion is accomplished 
through direct conversion of thermal energy to electric energy via thermoelectric 
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converters. The system's major components are the nuclear reactor, the primary 
coolant loop, the primary heat exchanger, the thermoelectric-electromagnetic (TE- 
EM) pump, and the thermoelectric converter. 

In the design presented by the GE Company [5], the lithium coolant exits the 
reactor s outlet plenum from six outlet lines. Each of the six outlet lines divides 
into two branches, and each of these twelve branches contains a TE-EM pump and a 
thermoelectric Power Conversion Assembly (PCA). After the lithium coolant exits the 
PC A, the twelve braches recombine into six reactor inlets. Each of the twelve primary 
heat transport branches is thermally coupled to a secondary coolant branch in the TE- 
EM pump and in the PCA where the thermoelectric elements are sandwiched between 
the primary and the secondary lithium coolant ducts. Each secondary lithium coolant 
branch is accompanied by a separate heat pipe radiator panel that radiates the waste 
thermal energy into space. Figure 1 is a drawing of the GE 100 kW e SP-100 system 
depicting the system’s key features [4]. 

II. 3 Statement of the Control Problem 

A major concern in the design of SNPSs is to ensure that these systems will 

operate safely and respond in a predictable and desired manner. Because of the 

design requirements of SNPSs and the stringent conditions under which they will 

have to operate, it is essential that the effects of an unanticipated system response 

be minimized. An unexpected system response may be due to degradation of the 
performance of the radiator heat pipes, loss of radiator heat pipes as a result of 
'ollisions with space debris, degradation of the performance of the thermoelectric 
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Figure 1. A Drawing of the General Electric SP-100 Thermoelectric Power System 
Design Depicting Its Key Features (Adapted from [6]). 


elements in the pumps and power conversion elements, fouling of the heat transport 
system, loss of a thermoelectric-electromagnetic (TE-EM) pump, sticking of the 
reactor control elements, and changes in the nuclear fuel characteristics due to fuel 
burnup, which might be reflected reflected in changes in the fuel temperature feedback 
coefficient . 
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Current designs of proposed static SNPSs require the use of banks of resistors 
(shunt regulators) to regulate the electric power delivered to the payload. Both 
thermoelectric and thermionic systems, employ a shunt regulator to reduce the power 
delivered to the load. In such systems, a failure in the shunt regulators could lead 
to reactor shutdown. Because of changing mission needs or deteriorating equipment, 
use of shunt regulators could also result in not delivering the expected electric power 
level, when the reactor operates below nominal thermal power. Furthermore, because 
the excessive load demand is rejected into space by radiation in the form of waste 
thermal energy, such an operation of the power system would increase the cooling 
requirements and hence the total mass of the spacecraft. In various studies by 
previous researchers, it has been demonstrated that a failure in the shunt resistors 
can not be tolerated [3], [4], [13]. Additionally, the results of studies investigating the 
load-following characteristics of static SNPSs without a shunt regulator show that 
although the nuclear reactor is always load-following, the integrated power system 
becomes nonload-following beyond a critical value of the load demand [4]. Therefore, 
a controller is needed in static SNPSs, to compensate for the possible reliability 
deficiencies resulting from the use of shunt regulators for the load- following of the 
power system. Such a controller could increase the overall system reliability, and 
it could improve the safe operation of both thermionic and thermoelectric SNPSs 
throughout their lifetime. 

11.4 Description of the Proposed Dual-Loop Controller 

A dual-loop control mechanism was initially proposed to eliminate the pos- 
sibility of a single-point failure in the shunt regulators of the static SNPSs. The 
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dUaUo °P con,roller ' whlch combines the advantages of an adaptive controller (inner 
loop). With those of an external controller (outer loop) is expected to demonstrate 
the load-following capability of a 300 kWe thermoelectric SNPS durmg its lifetime, 
by appropriately manipulating the neutronic, and thus the thermal power produced 
by the nuclear reactor. The controUed variable is the load voltage, regulated at 100 
V, whereas the manipulated variable is the core external reactivity, altered through 
the control drum motion. The dual-loop controller anticipates readings from sensors 
measuring the load voltage, the neutronic power, the primary coolant mass flow rate, 
the primary coolant temperature at the core inlet and exit, and the control drum 
position. Though in the final implementation of the external controller only the load 
voltage sensor reading is utilized. The controller generates a command for the new 


control drum position in order to meet the demanded electric power, while regulating 
the load voltage at the desired level. A block diagram of the dual-loop controller is 


depicted in Figure 2. 


II. 4.1 


Adaptive 

System 


Control Loop (Inner Loop) of the Static Space Nuclear Power 


The adaptive controller uses a model reference adaptive algorithm, which tends 
to minimize the difference between the output of a reference model and that of the 
actual system. The adaptive controller, which is activated only when the power 
system exhibits significant variations from nominal operating conditions, ensures that 
such variations do not hinder the performance of the external controller. The inner 
loop accomplishes this by causing the actual system response to match the expected 
response of the reference model, regardless of the complexity of the reference model 



A V*/ 



Figure 2. A Block Diagram of the Dual- Loop Controller 
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, . nUnt The design of this controller is fully described 

which embodies the physics of the plant. I he design 
in Metzger et al. [5]. 

„. 4 . 2 Nonlrnenr Extern*. Centro, Loop (Outer Loop, of .he S.atic Space Nuclenr 
Power System 

The nonlinear external controller forms the outer loop of the dual-loop con- 
ttoller. and is des.gned to operate when the power system responds as anticipated. 
The external controller is designed by gain-scheduling the linear controlle 
based on the reference mode, of the plant, and are obtained by using the LQG/LTR 
method. The outer loop, which is active a, all times is anticipated to exhibit Urn- 
jted robustness, and it is anticipated to be capable of overcoming system variations 
resulting from operation at different power levels. 

II. 5 Current Practice 

The voltage regulation problem in static SNPSs has not yet been solved. As ex- 
plained earlier, the current designs of the proposed static SNPSs using thermoelectric 

or thermionic energy conversion employ a shun, regulator (resistors) to reduce the 

is operating below its nominal thermal 


power 


delivered to the load when the reactor 


the load-following capabilities of 


power level. The results of the studies investigating 
i he TE power systems have demonstrated that : 

1 | A failure in the shunt regulator of the power system is a single p 

•2, The use of a shunt regulator for the load-following of the power system cannot 

be eliminated since the integrated power system operating without 
becomes nonload-following beyond a critical value of the load resistance ( 



18 


load resistance is equal to the internal resistance of the PCA), even though the nuclear 
reactor is always load-following, and, 

3 ) I'sing an adaptive control algorithm with selective state- variable weighting allows 
the power system to respond in a predictable, desirable and stable manner despite 
changes in the system's operating parameters. 

Even though reliability and vulnerability studies of static SNPSs have achieved 
a good understanding of the operating characteristics and the overall performance 
of the integrated system, they do not propose an effective control mechanism which 
could eliminate the possibility of a single-point failure resulting from failed shunt 
regulators. 
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CHAPTER III 

STATIC SPACE NUCLEAR POWER SYSTEM MODELING 

III. 1 Introduction 

The performance of the model-based controller proposed in this research de- 
pends particularly on the accuracy of the SNPS model used to desing it. Therefore, 
it is of great importance to obtain a model that embodies the physics of the SNPS un- 
der consideration. In this chapter, the description of an integrated, nonlinear model 
of the SNPS is presented along with a simplified, though nonlinear, reference model 
derived from it. The inputs, the outputs and the states of these models are described. 
The differences between the equations of the reference model and of the integrated 
SNPS model are pointed out. The linearization of the nonlinear reference model to 
obtain the linearized reference models about different equilibrium points is explained, 
and the range for which each linearized model is valid is determined. Then the lin- 
earized models which cover the entire operating regime (50% to 150% of the nominal 
-lectnc power level) of the power system are scheduled (fitted) as a function of the 
•dectric power produced by the system. An analysis of the linearized reference models 
and the scheduled reference model is presented and their validity is investigated via 
'tansient and steady-state response comparisons with the integrated, nonlinear SNPS 
model A FORTRAN computer program which incorporates the integrated, nonlinear 
SNPS model and the nonlinear reference model is used for these simulations. 
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III. 2 Description of the Integrated Nonlinear Space Nuclear Power System Model 

The integrated nonlinear SNPS model used in this research, was developed 
for analyzing integrated thermoelectric systems i4]. This coupled model is capable of 
predicting the performance of an integrated thermoelectric SNPS during both steady- 
state and transient operations, as well as predicting the system s load-following char- 
acteristics and the integration of the power system with the payload. Additionally, 
this power system model can be used to analyze the system response to various ex- 
ternal perturbations, the effect of system parameter changes on the transient system 
response, and the effectiveness of a power system controller. 

The major components of this system model are the nuclear reactor core, the 
primary coolant loop, the TE energy converter, the TE-EM pumps, the primary heat 
exchangers, and the main radiator. These components are listed in Table 1. A 
schematic diagram of the integrated, nonlinear SNPS model is depicted in Figure 3. 

The integrated nonlinear SNPS model has 16 states as shown in Table 2. 
The states of the model are coupled through a set of time varying, first order, 
nonlinear differential equations which represent the individual models that describe 
’lie integrated model. The equations that describe the neutronic power, the delayed 
neutron precursor concentrations, and the energy equations of the fuel, clad and 
’ht* reactor coolant are solved at each simulation time step using a stiff differential 
^nation solver. The remainder of the state equations are solved separately for each 
r 'nch at each time step also using a stiff differential equation solver. 
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Table 1. 


Major Components of the Integrated 


Nonlinear SNPS Model 


Component! 

Model 

Model Description 

Nuclear Reactor 

Neutron Kinetics and 
Reactivity Control Model 

* Point kinetics model with six delayed neutron groups. 

* Calculates the thermal power generated in the core due to fission 
for a step. ramp, or any arbitrary reactivity insertion function. 


Reactivity Feedback Model 

* Calculates reactivity feedback due to Doppler effect, fuel 
expansion, and coolant expansion. 


Reactor Thermal Model 

* A lumped parameter model. 

* Calculates the fuel and cladding average temperatures. 

* Calculates the core coolant average and exit temperatures. 

Reactor Hydraulic Model 

* Calculates the pressure drop in the reactor core region and in 
the annulus between the core baffle and the reactor vessel. 

Primary Coolant 
Loop 

TE-EM Pump Model 

• Couples the EM pump model to both the TE model and a 2-D 
thermal analysis model for the pump magnet 

• Calculates the pumping head and the overall pump efficiency. 

Piping Hydraulic Model 

* Calculates the pressure drop in the primary coolant loop pipings. 

Heat Exchanger Thermal 
Model 

* The energy balance model of the HX coolant channels is coupled 
to the main TE convener model. 

* Calculates the temperatures of the coolant and the wall in the HX 
coolant channels. 

* Calculates the total thermal energy extracted from die hear 
exchanger and the exit coolant temperature from the heat 
exchanger. 

Energy 

Conversion 

TE Model 

* A Transient 1-D finite element model. 

* Uses temperature dependent material properties. 

* Handles any combinations of boundary conditions. 

* Calculates the electric power output, conversion efficiency, 
and the waste heat rejection. 

Radiator 

i 

Waste Heat Rejection 
Model 

* Direct radiation model from the cold shoes of the TE generators 
into space. 

* Takes into account the effect of the sonic limit of the radiator s 
potassium heat pipes. 




22 



Figure 3. A Schematic Diagram of the Integrated, Nonlinear Space Nuclear Power 
System Model (Adapted from [6]). 

Ill 2 1 The Nuclear Reactor Model 

The neutron kinetics model solves the point kinetics equations using six delayed 
neutron groups. It calculates the time derivative of the reactor thermal power 
derated by fission, and the time derivative of the delayed neutron precursors. In 
model, the total reactivity input is the sum of the feedback reactivity caused by 
’^mperature variations (fuel expansion, Doppler, and coolant feedback effects), and 
external reactivity inserted into the core. 
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Table 2. States, Input, Output and Disturbance of the Integrated Nonlinear SNPS 
Model 


States 

Description 

Primary 

Neutronic power, MW t h 

Delayed-neutron precursor concentration I, MW t ^ 

Delayed-neutron precursor concentration II, MW ^ 

Delayed-neutron precursor concentration III, MW t h 

Delayed-neutron precursor concentration IV, M W t h 

Delayed-neutron precursor concentration V, MW ^ 

Delayed-neutron precursor concentration VI, MW t ^ 

Fuel temperature, K 

Average core coolant temperature, K 

Average primary branch coolant temperature, K 

Primary coolant mass flow rate, kg/s 

Secondary 

Thermoelectric temperature of the PCA, K 

Thermoelectric temperature of the TE-EM pump, K 

Secondary coolant temperature, K 

Secondary coolant mass flow rate, kg/s 

Radiator surface temperature, K 


Input 

External reactivity, A k/k 

Controlled Output 

Load voltage, V 

Disturbance 

Load demand, MW e 


The reactor core thermal model calculates the fuel and the average cladding 
temperatures, and the average core coolant and core exit temperatures using a 
lumped parameter approach. The model determines a centerline fuel temperature at 
each simulation time step, using the steady-state, one-dimensional heat conduction 
equation, assuming a flat radial power profile throughout the core. The average 
fuel temperature is calculated by averaging the clad temperature and the centerline 
temperature. The model calculates the core coolant exit and inlet temperatures as 
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the sum of the average core coolant temperature and the thermal energy added to 
or removed from the core coolant divided by the coolant mass flow rate through the 
core and the coolant specific heat, respectively. 

The core hydraulic model calculates the pressure drop in the reactor core as the 
sum of the pressure losses in the flow channels in the annulus between the reactor 
vessel and core baffle, and in the flow channels of the reactor core using the equivalent 
hydraulic diameter concept. 

III. 2. 2 The Primary Coolant Loop Model 

As shown in Figure 4, the primary coolant loop of the integrated nonlinear 
SNPS model consists of three flow paths and twelve heat exchangers with common 
inlet and exit headers. The coolant flows out of the core through three coolant flow 
paths, each equipped with a TE-EM pump. The coolant flowing out of the pumps is 
combined in the hot header and fed to the primary heat exchangers; the flow exiting 
the heat exchangers then combines in the cold header and returns to the reactor core. 
The TE-EM pumps are located in the hot leg of the primary coolant loop to maximize 
the TE converter’s efficiency and, consequently, the overall pump efficiency 4 . 

The TE-EM pump is a self-induced, direct current electromagnetic pump 
••quipped with a hiperco-27 magnet and a SiGe/GaP TE converter. Each pump has 
coolant ducts: the electric current provided by the pump s TE converter flows 
'■iirough the ducts in a perpendicular direction to both the magnetic field and the 
' "olam flow. Figure 5 depicts a cross sectional view of the TE-EM pump, and Figure 
’hows a line diagram describing the pump operation. The TE-EM pump model 
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Heat 



Figure 4. A Line Diagram of the Primary Coolant Loop (Adapted from [6]). 

assumes a constant coolant flow area and a uniform magnetic flux density through the 
coolant ducts. The model also assumes that the TE converter cold shoe temperature 
is equal to the pump radiator surface temperature. Thus, the model neglects the 
temperature drop between the evaporator and the condenser of the radiator’s heat 
pipes. Using these assumptions, the model determines the pressure rise across the 
pump duct, the electromagnetic (EM) pump efficiency, and the overall efficiency of 
'•he TE-EM pump. 

The piping hydraulic model calculates the pressure drop in the primary coolant 
°°P piping system. The model assumes the heat losses to be negligible, because 
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iTE-EM) Pump (Adapted from [6]). 
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Figure 6. A Schematic Representation of the TE-EM Pump Model (Adapted from 

[ 6 ])- 

'he primary coolant piping system is thermally insulated by insulation multifoils. 
( onsequently, the coolant temperature at the core exit can be assumed to equal that 
f the coolant inlet temperature to the TE-EM pumps. Additionally, the coolant inlet 
tnd outlet temperatures of the heat exchangers can also be considered equal to the 
'niperatures at the pump exit and at the reactor core inlet, respectively. 

The heat exchanger hydraulic model calculates the coolant flow rate in each 
'^Unt channel and the total pressure drop in the heat exchanger by ensuring that 
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the difference in friction pressure loss between any two consecutive channels equal 
the sum of the pressure losses in the hot and cold headers between the centers of the 
two channels. The model interpolates the friction factor as a function of Reynolds 
number. Such a function enables the integrated nonlinear SNPS model to respond to 
any changes in the reactor power and, consequently, to changes in the coolant flow 
from laminar to turbulent conditions, or vice versa. 

The heat exchanger model calculates the coolant and wall temperature distribu- 
tions in the coolant channels in the direction of the coolant flow. These temperatures 
are calculated through the coupling of the heat exchanger thermal model to the main 
TE converter model. 

111. 2. 3 The TE Converter Model 

The TE converter model allows the TE material properties such as the Seebeck 
coefficient, the thermal conductivity and the electrical resistivity to vary with tem- 
perature. The model calculates the electric voltage generated across the external load 
of the power system, the conversion efficiency of the TE generator and the electric 
power output of the system. 

111.3 Description of the Nonlinear Reference Space Nuclear Power System Model 

The reference model used in the design of the voltage controller is a simplified 
version of the integrated nonlinear SNPS model presented in the previous section, and 
it includes the physics that defines the system. The major differences between the 
nonlinear reference model and the integrated nonlinear SNPS model are the following 

[ 2 ]: 
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1. ) The nonlinear reference model has two delayed neutron precursor groups com- 

pared to the six groups of the integrated nonlinear SNPS model. The two-group 
delayed neutron precursor constants are obtained from the six-group delayed 
neutron precursor constants by using the method described by Skinner and 
Cohen [4j. 

2. ) The integrated nonlinear SNPS model provides for energy deposition in the 

nuclear fuel with a fuel clad, whereas the nonlinear reference model does not 
provide for a fuel clad. 

3. ) There are no momentum equations in the primary and secondary coolants loops 

of the nonlinear reference model. 

4. ) Values that can fluctuate in the integrated nonlinear SNPS model, such as 

the mass flow rate, remain constant in the nonlinear reference model for each 
operating power level. 

5. ) The heat transfer loop is modeled as a single lumped loop in the nonlinear 

reference model, and not as 12 separate branches as in the integrated nonlinear 
SNPS model. 

•5.) The secondary heat transfer loop and the radiator are lumped in one energy 
equation in the nonlinear reference model, while the integrated nonlinear SNPS 
model has a separate energy equation for each of the secondary coolant branches 
and each of the radiator panels. 

1 The nonlinear reference model does not have a separate set of energy equations 
for the TE-EM pumps as does the integrated nonlinear SNPS model, because 
the coolant flow rates are considered constant for each operating power level. 



Table 3. States, Input, Output and Disturbance of the Nonlinear Reference Model 


States 

Description 

Primary 

Neutronic power, MW t k 

Delayed- neutron precursor concentration I, MW^ 

Delayed-neutron precursor concentration II, MW t h 

Fuel temperature. K 

Coolant core exit temperature, K 

Heat exchanger exit temperature, K 

Coolant core inlet temperature, K 

Secondary 

Thermoelectric temperature, K 

Secondary coolant temperature, K 


Input 

External reactivity, &k/k 

Controlled Output 

Load voltage, V 

Disturbance 

Load demand, MW e 


The nonlinear reference model has nine states, one input, four measured outputs 
and one disturbance, as shown in Table 3, and it is described by the following nine, 
first-order differential equations: 

1.) The reactor kinetics equations : 

iMl = (i) 


where, 

P n {t) : the reactor neutron power (the neutronic power), 

C^\t) : the delayed neutron precursor concentration for the i tk group, 

p t (t) : the total core reactivity, 

flt : the toted delayed neutron fraction, 

A : the neutron generation time, and, 
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A'*' : the delayed neutron precursor constant for the group, 
where the delayed neutron precursor concentrations are determined from 

<fC (1) (0 = $^_p n{t) _ A (1) c {1] (t), (2) 

dt A 

and 

= —P n (t) - A (!I C (II «), (3) 

dt A 

where, 

: the delayed neutron fraction for the group, 
and, where the total core reactivity used in the reactor kinetics equations is given by 

pt(t) = Pez{t) + Qpin [ f7^To) ] + Q /«i[ r /(0 ~ T A l = °)1 + ^[^(O - T cr (t = 0)], 

(4) 

where, 

p ex {t) : the external reactivity, 
q £f ; the Doppler reactivity feedback coefficient, 

Tf(t ) : the average fuel temperature, 

a/ ue i : the fuel temperature reactivity feedback coefficient. 

a CT : the coolant temperature reactivity feedback coefficient, 

T cr (t ) : the average core coolant temperature. 

- 1 The nuclear fuel energy equation : 

dTfjt) Pn{t] {Tfjt) — O.oTcroU) - 0-5T tn (Q) { . 

dt M fC p f T f 


where. 
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A If : the mass of the nuclear fuel. 

c p f . the specific heat of the nuclear fuel. 

T C ro{t) : the core coolant exit temperature, 

T in (t) : the core coolant inlet temperature, and where tj is defined in the 
following paragraphs. 

3.) The cold-leg energy equation 

dT m (t) = (T hxo jt)-TM) (6) 

dt T( n 


where, 

Thxo(t) '■ ^e heat exchanger coolant exit temperature, and where r tn is defined 
in the following paragraphs. 

4.) The reactor coolant energy equation 


dT C ro(t) (^/(O — 0-5 T C ro(t ) 0.5 Ti n (t)) mtot 


dt 


Mr 


{T C ro(t) - Tin(t)), (7) 


where, 

M c T : the core coolant mass, 

m tot : the total core coolant mass flow rate, and where, r cr is defined in the 
following paragraphs. 

V i The primary loop coolant energy equation : 

dTk xo (t) mtot rp ,,, T ,.u (0.5 T cro (t) — 0.5 Th xo {t) — Tt e (t )) 

-r = -rr(T ero [t) - T hxo (t)) , (») 

dt M c T cp 

« here. 


Me the coolant mass in the primary heat transfer branches. 
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T„w : .he thermoelectric temperature, and where r cp , s defined tn the following 
paragraphs. 

6. ) The PCA thermoelectric energy equation : 

kZ T,kz MteCpte 

(9) 

where, 

Jt °t ( 0 ; current drawn by the load, 

^loadit) ■ the voltage drop across the load, 

Mtc . the mass of the thermoelectrics, 

: the specific hea, of the PCA .hermoe.ec, ties, and where r kt and r, k , are 
defined in the following paragraphs. 

7 ) The secondary coolant energy equation : 


_ (Tte(t) - T,(t)) (T,(t) - T sxnk ) 

dt -r I ! 


where, 


^(f) - secondary coolant temperature, 

T ‘"‘ t = ' ff ' C,iV5 smk l ' mperature - »l>ere, r„ and r. are defined in the 
following paragraphs. 

The r’s in the energy equations (Eqs. (5) through (10)) are „me constants that 

are .-..tally calculated for each operat.ng power level, from the phystcal propert.es 

of the system such tha, the equat.ons satisfy steady-state condit.ons. The r’s are 
defined as follows: 


T f = (T f (t) - 0.5 T cro (t) - 0.5r, n (t))— /Cp/ . 

1 Pn(t) ’ 


(11) 
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where 


T rr — 


{T f (t)-Q.hT C ro(t) -0.5T tn (f))M c 
(T C ro(t) - T in (t))rh tot 


( 12 ) 


T cp — 


(O. bjcrojt) + 0.5 T hxo (t) - r (e U))A/ c 

( T C ro{ t ) — Tl LX0 (t))TTlt o t 


(13) 


T hx — 


T cpM te Cp te 

M c c p c 


14) 


T skx ~ 


(T te (t)-T,(t)) 


£(0.5 T„ 0 (t) + 0.5 T hxo (t) - Tu(t)) - 


do) 


t, t = 


^~sh.x M,C 


pa 


MteCp te 


(16) 


and 


T r = 


Ta riUt) - T , ink (t)) 
{Tte{t) ~ T,(t)) ’ 


(17) 


Cpc : the primary coolant specific heat, 

M s : the coolant mass in the secondary coolant branches. 
c p} : the secondary coolant specific heat. 

Equations (1) through (10) can be placed in the following state-space variable 


X(f ) = f (x(t), u(t), u;(t )), 


whereas the output equation can be defined by 


(18) 


y (t) = g (x(t), u(t) , iv(t)), 


(19) 
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where f(.) and g(.) are vectors of nonlinear functions of the system states, input, and 
disturbance. The state vector x(t) is defined by 


x(<) 


Pn(t) C' l Ht) CW(t) Tf(t) 




T’cro(0 Tfizo(0 T in (t) Tte{t ) T s (t) , 


( 20 ) 


and, 

u(t) is the system input (external reactivity, p ex (t) ). 

w(t) is the system disturbance (the electric load demand, Itot{t) L/oad(0 )• and, 
y(t) is the measured system output (5 x 1). 

The measurable outputs of the system are as follows : 

1. ) The neutronic power, P n {t), 

2. ) The core coolant exit temperature, Tcro(t), 

3. ) The core coolant inlet temperature, Tin(t), 

4. ) The total core coolant mass flow rate, m(t), and, 

3.) The load voltage, Vi oa d{t)- 

Equations (18) and (19) are solved at each time step during the simulation 
using a stiff differential equation solver, as it is done when simulating the integrated 
uonlmear SNPS model. 


4 Linearized Reference Models 

• 11 t l Development of the Linearized Reference Model 

The linearized reference models are obtained for each operating power level 
• lna lvtically and numerically linearizing the nonlinear reference model. Analytical 
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linearization is performed on the nonlinear reference model of the SNPS using a Taylor 
series expansion about the stead-state values of the states, input and disturbance 
(x 0 , u 0 , u; 0 ). These steady-state values are obtained from operating the integrated 
nonlinear SNPS model at various power levels. Applying the Taylor's expansion to 
Equation (18), and neglecting higher order terms, the following is obtained: 

df 

x 0 (0 + $x(0 =f [xe(i)ti*o(t)'WD(t)] + l(x 0 .u 0 .wo)* x (0 

+ (^) i(*o,»o,-o)MO + ) l(x 0 ,tto.u>o)^ u/ (*) ^1) 

where 5(.) is the perturbation from the equilibrium point and where all the derivatives 
are calculated at each operating power level. Because the system is at equilibrium at 
(Xo, uo, u>o), Equation (21) can be rewritten as : 

+ (£) !(«.»«)' M«) + (£) < 22 > 

As a result of the hard nonlinearities present, the output Equation, (19), is 
numerically linearized by perturbing the system states, input variable and disturbance 
enable about an equilibrium point (steady-state conditions corresponding to a 
particular power level) and the rates of change of the system states and the output 
enables are determined. Expanding Equation (19) in a Taylor series about the 
quihbnum values of the states, input and disturbance, the following is obtained: 



(23) 


* l (xo.uo.u. 0 )^ u ( < ( * ■{x OW ’0) 6w ( t ) 

where SI.) is the perturbation from the equilibrium point, and where all of the 
derivatives are evaluated at (x 0 ,u 0 .w 0 ). In Equation (23). the output refers to the 
load voltage alone, because only one output can be independently controlled by the 
single input control. 

From Equations (21) and (23), the following linearized reference system model 
is obtained which is valid about the equilibrium point (x = xo, u= uq, u>=it/o) : 


6x{t) = A6x(t) + B6u(t) + L6w(t) (24) 

<5 y(t) = C6x(t) + D<5u(f) + F£u/(<) (25) 

where £x(t), 6u(t), 6iv(t), and Sy{t) are the vectors perturbed about their equilibrium 
values. Figure 7 depicts a block diagram of the linearized reference model. 

The system matrices are given by the following expressions: 


A=( s> 

'■ ( Xo.UQ.U’o) 

, State matrix (9x9) 

(26) 

»-£> 

( Xo.uq.u;o) 

, State Input vector (9x1) 

(27) 


I Xo.Uo.U’o) 

. State Disturbance vector (9x1) 

(28) 

C = ( — ) 

dx 

( Xo.U(),U*o ! 

, Output vector (1x9) 

(29) 

d = ( f> 

au 

{ i Xo.^0« u, 0* 

, Feedforward scalar term 

(30) 

OlV 

'( Xo,uo*wo} 

, Output disturbance scalar term 

(31) 



&c(t) : State Vector 
8u(t) : Input Vector 
8y(t) : Output Vector 
5w(t) : Disturbance Vector 


?ure 7. Block Diagram of the Linearized Reference Model 
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The elements of the linearized reference model are as follows : 





4i2 = Ai, 


^13 = 

A\i = 
A 21 = 
-422 = 
A31 = 
•432 = 
•441 = 

444 = 

445 = 
-447 = 

4$4 = 
455 
457 
•466 

467 

468 
476 

4 77 


h, 


<*ez 

A ’ 

§1 
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-Ai, 


02 

T’ 


-^ 2 , 

Pn(t=0) 
^t c v! ' 

-1 
r f ' 

05 
r t ’ 

0.5 
T / ’ 

1 

Ter' 


—0.5 mtot 


r C r 


_ ™tot 

m 7 ’ 


—0*5 , rhtot 
T C r + Me » 

mtot 0.5 

MT r cp ’ 

0.5 


= lf- 


T C p’ 


__ 1 


__ 1 


> cp 5 


- <Li 

r fcx’ 





^86 = 


0.5 


r hz ’ 

^88= r 1 -, 


T hx 


' thx 


A 89 = ^7’ 

^ 98 = 



J_ 

TV 1 


#11 = x 


181 - 

where, A is a (9 x 9) matrix, and B and L are (9x1) column vectors. 

The numerical values of the elements of A, B. C, L, and F for various 
equilibrium points are given in Appendix A. All other elements not listed above 


or in Appendix A are identically zero. 
III. 4. 2 Linearization Error Analysis 


The linearized reference model obtained by linearizing the nonlinear reference 
model at a particular operating point is valid only for a small perturbation around this 
operating point. Therefore, a linearization error analysis is required to determine the 
magnitute of the perturbation about the equilibrium point for which the linearized 
model is valid. This range of validity is determined by comparing the state and 
output (load voltage) responses of the linearized reference model with those of the 


integrated nonlinear SNPS model for various perturbations in the external reactivity. 


corresponding to equivalent changes in the load demand. The percentage relative 
error (e) in the states of the linearized reference model is defined as : 


£ 


(fo'niU) 

max 

t £ tmor 



i = l 9 


(32) 
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where 6x l nl (t ) and 6x'i(t) are the state perturbations of the nonlinear integrated SNPS 
model and linearized reference models from the equilibrium point, respectively, and 
x l ni(t) is the absolute state of the integrated nonlinear SNPS model. In this analysis, 
tmax is assumed to be 400 seconds since this covers the majority of the transient and 
steady state responses of the system. Step perturbations of the external reactivity 
from equilibrium, both in the positive and negative directions, are applied until the 
maximum percentage relative error, e. of any of the states and the load voltage reaches 
about 10%. The maximum power level perturbation (positive or negative) for which 

this value of e is reached is chosen as the validity range for the linearized reference 
model. 

To obtain the models that cover the entire operating regime of the SNPS (from 
50% to 150% of the full power), a linearized reference model of the SNPS is obtained 
at 100% of full power level. A linearization analysis is performed for this model and 
its validity range in the positive and negative directions is determined. The next 
operating point about which a linearized reference model is obtained is determined 
by considering the validity range of the previous model, and the extent to which two 
consecutive models will overlap. Because the validity range of the next model is not 
known apnori, this process is iterative in nature. The iterations are terminated when 
satisfactory overlapping is achieved. The linearization process is repeated about each 
selected operating point both above and below the full power level using the above 
procedure until the entire operating regime of the power system is covered. 
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Table 4. Validity Ranges of the Linearized Reference Models for Positive and Neg- 
ative Step Perturbations in the Electric Power Output 


Equilibrium 
Power Level 

Maximum 

Positive 

Perturbation 

Maximum 

Negative 

Perturbation 

50% 

5% 

5% 

55% 

5% 

5% 

60% 

5% 

5% 

65% 

8% 

6% 

70% 

10% 

8% 

75% 

15% 

12% 

80% 

18% 

15% 

85% 

19% 

15% 

90% 

23% 

20% 

95% 

29% 

26% 

100% 

38% 

24% 

105% 

38% 

23% 

110% 

37% 

21% 

115%* 

35% 

19% 

120%* 

30% 

17% 

125%* 

25% 

17% 

130%* 

20% 

15% 

135%* 

15% 

10% 

140%* 

10% 

6% 

145%* 

5% 

5% 

150%* 

5% 

5% 


- shows the models are valid at and beyond 150% 


Table 4 depicts the validity ranges for the linearized reference models at different 
equilibrium points. It is observed that the validity ranges of the linearized reference 
models are quite large and the amount of the perturbations for models above 100% of 
full power are larger than those for the models designed for operating levels below full 
power. This can be explained by the fact that when the reactor operates above full 
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power, the responses of the system states and the output are varying less from the 
nominal power level than when the system operates below the nominal power level, 
due to a higher efficiency of the TE-EM pumps at power levels above the full power 
level. Therefore a smaller change in the external reactivity is required to bring the 
system at full power to a higher equilibrium point than to bring it to an operating 

level below the full power level. 

HI. 5 Validation of Linearized Reference Models 

In this section, a comparison of the state and output (load voltage) responses 
of various linearized reference models and the integrated nonlinear SNPS model to 
certain step changes in the external reactivity are presented. Figures 8 through 15 
depict the open-loop responses of the load voltage and of the states of both the 110% 
of full power linearized reference model and the integrated nonlinear SNPS model to 
a step increase in the external reactivity, corresponding to a step load level increase 

from 100% to 110% of full power. 

Similarly, Figures 16 through 23 depict the open-loop responses of the load 
voltage and of the states of both the 80% of full power linearized reference model and 
the integrated nonlinear SNPS model to a step load level decrease from 100% to 80% 

of full power. 

Figures 24 and 25 show the responses of the neutronic power and the load voltage 
of the 110% of full power linearized reference model and the integrated nonlinear 
SNPS model to a step increase in the external reactivity, corresponding to a 35% of 
full power increase in the load level of the system, respectively. It is observed that 
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nnl 
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Figure 8. Load Voltage Responses for a Step Change in External Reactivity Corre 
sponding to a Load Change from 100% to 110% of Full Power. 
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Figure 10. Fuel Temperature Responses for a Step Change in Reactivity Correspond- 
ing to a Load Change from 100% to 110% of Full Power. 
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Figure 11. 


Core Coolant Exit Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 110% of Full 
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Time (sec) 

Figure 12. Heat Exchanger Coolant Exit Temperature Responses for a Step Change 
in External Reactivity Corresponding to a Load Change from 100% to 
110% of Full Power. 
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Figure 13. Core Coolant Inlet Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 110% of Full 
Power. 
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Figure 14. Thermoelectric Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 110% of Full 
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Figure 16. Load Voltage Responses for a Step Change in External Reactivity Cor 
responding to a Load Change from 100% to 80% of Full Power. 
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Figure 18. Fuel Temperature Responses for a Step Change in External Reactivity 
Corresponding to a Load Change from 100% to 80% of Full Power. 
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Figure 23. Secondary Coolant Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 80% of Full 
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the output and state responses of the linearized reference model match the integrated 
nonlinear SNPS model responses quite closely. 

III. 6 Scheduling of Linearized Reference Models 

Scheduling refers to the fitting of the linearized reference models to a variable 
which is indicative of the system changing operating conditions. In this case, the 
variable chosen is the electrical power (P e ) generated by the system. The scheduled 
plant dynamics can be represented as: 

6x(t) = A(P e )6x(t) + B(P e )M0 + L(P e )M0> (33) 

Sy(t)= C(P t )6x(t) + D(P e )6u(t) + F(P e )6w(t), (34) 

where Sy(t ), <5x(<), 8u(t), 8w(t) are the output, states, input and disturbance vectors 
of the scheduled reference model, respectively. A, B, F, C, D, and L are the 
scheduled model system matrices with elements that axe functions of P e . It is 
anticipated that the scheduled model reference will be more accurate than any single 
linearized reference model in the entire operating range of the SNPS. 

III. 7 Evaluation of the Scheduled Reference Model 

Figures 26 through 33 show the responses of the load voltage and of the states of 
the 80% linearized reference model, the scheduled reference model and the integrated 
nonlinear SNPS model to an external step reactivity insertion of -9.51, cents which 
corresponds to a load level change from 100% to 80% of full power. Similarly, 

Figures 34 through 41 compare the responses of the load voltage and the states of 
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Figure 24. Neutronic Power Responses for a Step Change in External Reactivity 
Corresponding to a Load Change from 100% to 135% of Full Power. 
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Figure 25. Load Voltage Responses for a Step Change in External Reactivity Cor 
responding to a Load Change from 100% to 135% of Full Power. 
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Figure 26. Load Voltage Responses for a Step Change in External Reactivity Cor 
responding to a Load Change from 100% to 80% of Full Power. 
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Figure 27. Neutronic Power Responses for a Step Change in External Reactivity 
Corresponding to a Load Change from 100% to 80% of Full Power. 
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Figure 29. Core Coolant Exit Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 80% of Full 
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Figure 30. Heat Exchanger Coolant Exit Temperature Responses for a Step Change 
in External Reactivity Corresponding to a Load Change from 100% to 
80% of Full Power. 
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Figure 31. Core Coolant Inlet Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 80% of Full 
Power. 
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the 150% of full power linearized reference model, the scheduled reference model and 
the integrated nonlinear SNPS model to +23.73 cents of reactivity insertion, which 
corresponds to a load level change from 100% to 150% of full power. 

The response of the secondary coolant temperature of the linearized reference 
model shows variations from that of the scheduled reference model; the scheduled 
reference model improves the steady state behaviour of this state in the expense of 
causing a non-minimum phase behaviour in the transient response of the state, even 
though the integrated nonlinear SNPS model is minimum phase. The deviation of 
the secondary coolant temperature response of the scheduled reference model from 
the integrated nonlinear SNPS model response has been determined to be less than 
0.01%, therefore this non-minimum phase behaviour can be attributed to modeling 
error. By examining the aforementioned figures, it is observed that the scheduled 
reference model matches the integrated nonlinear SNPS model states better than a 
linearized reference model. 
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Figure 34. Load Voltage Responses for a Step Change in External Reactivity Cor- 
responding to a Load Change from 100% to 150% of Full Power. 
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Figure 35. Neutronic Power Responses for a Step Change in Negative Reactivity 
Corresponding to a Load Change from 100% to 150% of Full Power. 
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Figure 36. Fuel Temperature Responses for a Step Change in External Reactivity 
Corresponding to a Load Change from 100% to 150% of Full Power. 
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Figure 37. Core Coolant Exit Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 150% of Full 
Power. 
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Figure 38. Heat Exchanger Coolant Exit Temperature Responses for a Step Change 

in External Reactivity Corresponding to a Load Change from 100% to 
150% of Full Power. 
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Figure 39. Core Coolant Inlet Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 150% of Full 
Power. 
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Figure 40. Thermoelectric Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change from 100% to 150% of Full 
Power. 
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Figure 41. Secondary Coolant Temperature Responses for a Step Change in External 
Reactivity Corresponding to a Load Change Rom 100% to 150% of Full 
Power. 
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CHAPTER IV 

LINEAR CONTROL SYSTEM 
DESIGN AND ANALYSIS TOOLS 


IV. 1 Introduction 

The linearized reference models of the SNPS were found to be open-loop stable, 
yet load voltage responses of the linearized referencemodels have very short rise times. 
Thus, it is desired to design a compensator that : 

1. ) Slows down the load voltage response of the system, allowing more time for the 

insertion or removal of reactivity by the control drums. 

2. ) Achieves desirable performance such as command-following, disturbance rejec- 

tion, sensor/actuator noise insensitivity at high frequencies, and, 

3. ) Exhibits stability robustness to system parameter variations, unmodeled dy- 

namics, etc. . 

In this chapter, some of the mathematical preliminaries, such as the definition 
and some properties of singular values, are discussed. Also, some loop shaping 
concepts are described in relation to the singular values of system transfer function 
matrices. The multivariable linear control system design method, LQG/LTR, is 
described based on a Model-Based Compensator (MBC) structure. Some of the 
command following, disturbance rejection and robustness properties of the LQG/LTR 


compensators are also presented. 
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At this stage it is pertinent to remark that even though the LQG/LTR tech- 
nique is for MIMO systems, it is used to design Single Input-Single Output (SISO) 
compensators for voltage regulation of the SNPS system, primarily for the following 
reasons: 

1. ) There are several desirable guaranteed properties of the LQG/LTR compen- 

sators, such as gain and phase margins, and stability robustness, and, 

2. ) The design of an LQG/LTR compensator is done in a systematic manner, 

allowing the designer to address issues such as system bandwidth and other 
frequency domain specifications in a more direct manner than what is possible 
using conventional SISO control design methods, while performing the design 
using time-domain tools. 

IV.2 Mathematical Preliminaries 

Some of the key results and design considerations of the LQG/LTR technique 
are stated in terms of singular values of the various system matrices. The singular 
values of a complex n x m matrix A, denoted as <ri(A), are the k largest nonnegative 
square roots of the eigenvalues of (A H A), where k = min(n,m). That is : 

<r;(A) = y/a,(A h A), i = 1,2, . . . , fc, (35) 

where A® is the complex conjugate transpose of A. One representation of the matrix 
A, known as the singular value decomposition (SVD), is given by: 

n 

A = U £ V H = <Ti{ A) u; vf, 

1=1 


( 36 ) 
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where. 




u 2 


l>2 


j; U H U = I, 
V H V = I, 


V = 

E = diag (<ri,<T 2 , 


(37) 

(38) 

(39) 


where the columns of the matrices V and U are the eigenvectors of A® A and A A**, 
respectively. The maximum and minimum singular values may also be defined in 
terms of the matrix spectral norm ||-|| of the matrix, as follows [17] : 

= IST = l|A|1 ” (40) 

and 

tfmtn = min = ||A -1 || , (41) 

x#o ||x || 2 11 1,2 v ’ 

Also, it is true that 

Vmxn = j.J| A 1 H 2 ’ ifdetA^O (42) 

l 0, if det A = 0. 

A useful property of singular values is the following inequality: 

<7m«n(A) > <Tmar(B) =* A® A > B® B. (43) 


Additionally, the triangle inequality holds good for the maximum singular value. That 
is, 


Gmax (A + B) < ^maz(A) ■+■ 0rruiz(B). 


( 44 ) 
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Also, if the inverse of the matrix A exists, then from Equations (40) and (42) the 
following is obtained: 


^mtn( A ) 


1 

^’mar(A. *) 


(45) 


IV. 3 Loop Shaping Concepts 

Generally speaking, loop shaping refers to the design of control systems by 
influencing the feedback characteristics to correspond to a specific form which is 
known to result in systems with desired properties, such as closed loop stability, 
command following, disturbance rejection, and insensitivity to sensor and/or actuator 
noise. In the framework of SISO systems, the quantification of the characteristics of 
the Transfer Function Matrices (TFMs) such that the resulting closed-loop system has 
desired properties is fairly straightforward [18], [19]. The block diagram of a general 
MIMO system is shown in Figure 42. From Figure 42 it can be concluded that the 
closed-loop system output is given by the following expression : 

y(s) = [I + G{s) K(j)f \ G(j) K(«) r(«) + [I + G(a) K(a)] -1 d 0 (s) 

+ [I + G(s) K ( a)]” 1 G(s) di(*) - [I + G(s) K(s)] -1 G{s) K(a) n(s), (46) 

where, 

r(s) : the reference input vector (m x 1), 
dj( 5 ) : the input disturbance vector (m x 1), 
d 0 (s) : the output disturbance vector (m x 1), 
y(s) : the output vector (m x 1), 
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The general design requirements for a MIMO feedback loop are: a) Closed-loop 
stability, b) Good command following, c) Good disturbance rejection, d) Sensitivity 
reduction to unmodeled dynamics and e) Sensor noise insensitivity. These design 
requirements can be related to the system TFMs as follows: 

1). Command Following 

Command following is the ability of the system output y(t) to follow the 
reference input r(£) closely. In order to achieve good command following it is required 
that: 

y(s) ~ r(s) for all s € fi rt (47) 

where fi r is the frequency range of expected references. From Equation (46), the 
following condition for good command following is obtained : 

[I + G(s) K(s)]" 1 G(s) K (s) ~ I, s € (48) 

which can be accomplished if the return difference TFM has the following property: 

[I + G(s)K(s)] is “large” for s € (49) 

which implies that the forward loop TFM 

G(j) K(s) is “large” for s 6 H r . (50) 

Thus, for good command following, the transfer function matrix G(s)K(s) must be 
“large” for the expected range of reference input signal frequencies, provided that the 
closed-loop system is stable. 
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Disturbances 



Figure 42. Schematic Diagram of a General Multiple Input-Multiple Output 
(MIMO) System. 


n(s) : the sensor noise vector (m x 1), 

and G(s), K(s) are the (m x m) plant and compensator TFMs, respectively. From 
Equation (46), it can be inferred that: 

a. ) [I + G(s) K(s)] -1 G(s) K(s) relates the reference input to the output, 

b. ) [I -I- G(s) K($)] -1 G(s) relates the input disturbances to the output, 

c. ) [i+g^koof 1 di(«) relates the output disturbances to the output, and, 

d. ) [I + G(s) K(s)] -1 G(s) K(s) relates the sensor noise to the output. 
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2.) Disturbance Rejection 

Disturbance rejection is the insensitivity of the closed-loop system to input 
disturbance signals, dj(s), d 0 (s). In view of the discussion on good command following 
properties of the system, and from Equation 46, it can be observed that the output 
y(s) is influenced by disturbance signals which could force it to deviate from the 
reference, r(s). In a manner similar to the development above, it can be shown that 
disturbance rejection is accomplished if: 

[I + G(s) K(s)] is “large” for s e s e Q do , (51) 

where Q di and Q do are the frequency ranges of the expected disturbances, or, equiva- 
lently 

G(s) K(s) is “large” for s € s 6 Q do . (52) 

Equations (51) and (52) require that the forward loop TFM must be large in the 
frequency range where the disturbances have most of their energy. 

3.) Sensitivity Reduction to Unmodeled Dynamics 

Sensitivity reduction of a control system refers to its ability to attenuate the 
effects of deviations of the plant from its nominal values as well as the effects of 

unmodeled dynamics. Consider the nominal output y*( a ) of a MIMO system to a 
reference input r(s) : 


y'(«) = [I + G'( 3 )K(j)1 _i G»K(,)rM, 


(53) 
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where G‘(s) refers to the nominal open-loop system. Let the actual open-loop system 
be given by: 


G(s) = G*(s) + £G(s), 

and the actual output by: 

y(') = y» + Sy(s). 

Now, from the above equations, 


(54) 


(55) 


y-( S ) + Sy(s) = 11+ [G» + «G( S )] K( S )]-‘ [G-(») + «G(,)] K( 5 ) r(,). (56) 


From the above equations, it can be shown that: 


W*) = ( J + G (*) K(^)]- 1 SG(s) G'(s)- 1 y *(,). (57) 

From Equation (57) it is observed that if the loop TFM, G(s)K(s), is chosen to be 

large (and provided that the system is stable) the effects of SG(s) on 6y(s) can be 
reduced to a great extent. 

4.) Sensor Noise Insensitivity 

From Equation (46) and from the above arguments it can be easily seen that the 
system will be insensitive to sensor noise by ensuring that the loop TFM, G(s)K(s), 

is small in the frequency range where most of the sensor and actuator noise exists i.e. 
at high frequencies. 


IV.4 The Role of Singular Values in Loop Shaping 

Most of the loop shaping concepts, as discussed above, can be reduced to 
specifying the magnitude of the forward loop TFM, G( j) K (*), or the return difference 
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TFM, [I + G(s) K(s)], in the specific frequency ranges of interest. It is in this context 
that the concept of singular values plays a role. The relations between the maximum 
and minimum matrix spectral norm and singular values (<r) are as follows : 


max 

u(a)^0 


jgMfllla _ , n , » 

— 0rna;c( G($)), 


and, 


min 


l»M 


|G(,) «(,)!]*_ 

— 0Vntn(G(s)). 


(58) 


(59) 


u( s )#o ||u(a)j| 2 

Thus, the transfer function matrix G(«) is said to be “large”, if <r mtn (G(s)) is large. 
Similarly, G(s) is said to be “small”, if <r Tnar (G(s)) is small. 

Therefore, MIMO loop shaping concepts can be restated in terms of the sin- 
gular values of the loop TFM and the return difference TFM That is for good com- 
mand following, disturbance rejection, sensitivity reduction, and noise insensitivity 
0mtn(G(s)K(s)) should be large in the low frequency region, and < 7 mar(G(s) K(s)) 
should be small in the high frequency range. 


In addition to the above relations in which the system performance specifications 
are stated in terms of singular values of the appropriate TFMs, it can be shown that 
the singular values of the appropriate transfer function matrices are good measure of 
robustness to modeling errors [20]. This is explained as follows: if E(s) is the error 
between the actual system denoted as G(i) and a model of the system, G(s), then 
there are relations in terms of singular values of the system that guarantee closed- 
loop stability of the actual system when the compensator is designed using the model 
G(s). As the error E(s) is not known exactly, a conservative bound is used. 
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Define a matnx L(s), as a measure of uncertainty of the model, G(s), of the 
system. Therefore, the relation between the actual plant, G(s), and the model, G(s), 
is: 

G(^) = L(s)G(s). (60) 

Define E(s) such that the following relation is satisfied: 

L(s) = E(s) + I. (61) 

From Equations (60) and (61) the following expression for E(s) can be derived: 

E(j>) = [G($) - G(s)] G _1 (s). (62) 

It has been shown by Lehtomaki [20] that a compensator K(s) which stabilizes the 
closed-loop based on the model G, will also stabilize the closed-loop with of actual 
plant G, if the following relation is true: 

~ I] < a ’min [I + (K(joi) G(jfai)) -1 ] . (63) 

It can be shown that at high frequencies (when most of the modeling errors are 
apparent), if |L(;u,)| » 1 and | (G(ju,) K(ju)) | » 1, then Equation (63) becomes 

o-mai [L(ju,)] < a mxn [(K(>w) G(>w)) -1 ] = — - 1 (64) 

&maz [K(ju;) G(ju;)] 

Thus, by using Equation (64) and the high frequency noise rejection specification 
that <r mttI (G(s)K(s)) is “small” in the high frequency range, the high frequency 
uncertainly requirements for the system model can be defined. 

Thus, the above discussion summarizes the relation between the MIMO loop 
shaping concepts and the singular values of the system TFMs. 
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IV. 5 Model- Based Compensation 

Dynamic compensators are required for plants which must exhibit different 
behavior at high and low frequencies. Model-Based Compensators (MBCs) are a 
special class of dynamic compensators used in feedback control when only part of the 
state is measured. Use of MBCs allows the indirect estimation of Jill the unobservable 
states for use in the feedback control. A block diagram of an MBC in a closed-loop 


is shown in Figure 43. 

The state-space representation of the plant is given by: 

x(() = Ax(t) + Bu(t) + Lw(t) ; x(t) € » n ,u(t) € R m , w(t) € (65) 

y(f) = Cx(t) ; y(t)€» m , (66) 

whereas the state-space representation of the MBC is given by: 

z(t) = Az(l) + Bu(t) + H[-e(0-Cz(t)] ;z(t) € » n ,e(t) 6 » m , (67) 

u(i) = — Gz(t), (68) 

or, by, 

z(t) = [A-BG-HC] z(t)-He(t), (69) 

u(0 = -Gz(t), (70) 


where z(f) represents the MBC state vector. The feedback interconnection relation 
is given by: 


e(() = r(e) - y(<) ; r (t) 6 S", 


(71) 
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where r(t) represents the reference vector. The state vector for the closed-loop system 
is defined as: 

*cl(0 = MO *(0]' ; Xcl(t) € tf 2n , (72) 

where 

x(0 = x(t)-z{t). (73) 


The state-space representation of the closed-loop system can now be expressed as : 

MO 

MO 

A c j 


A - BG 

0 


-BG 

A - HC 


x(0 

x(0 


+ 


0 

-H 


r(0 + 


0 

-L 


w(0, (74) 


y(0 = [c 



(75) 


Additionally, it can be shown that: 


det(A I - A d ) = det(A I - A + B G) . det(AI - A + H C). (76) 

From the above representation of the MBC, H and G are the design matrices which 
must be chosen such that the following conditions for nominal stability of the closed- 
loop system are satisfied: 


ReA;(A-BG)<0 ;i = l,2,...,n, (77) 

Re A, (A — H C) < 0 ;i = 1,2,... ,n. (78) 

In addition to ensuring nominal plant stability, the MBC yields systems with special 
properties; the design matrix G controls n poles while the design matrix H controls 
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the other n poles. By virtue of this separability of the closed-loop poles, the poles 
of (A - BG), can be placed using full-state feedback designs such as the Linear 
Quadratic Regulator (LQR) or even pole placement techniques [21]. This is known 

as the control loop design and its state-space representation (the full-state regulator 
problem) is: 


x(t) = Ax(t) + Bu (t), 

u(t) = -Gx(t). 

The state-space representation of the closed-loop system is given by 


(79) 

(80) 


x(t) = [A - B G] x(t). (81) 

The other half of the closed-loop poles are the poles of (A - HC), and they can 
be placed using the Kalman Filter (KF) algorithm to guarantee nominal stability. 
This is known as the filter loop design and its state-space representation (an output 
feedback problem with no external control) is: 


x(0 = Ax(t) + v(f), 


(82) 


y(t) = Cx(<), 


(83) 


Wh " e V(<) repr “ enls ,he '“novations term. This term is related to the sensed ontput 
as follows: 


v(t)= -Hy(<) = -HCx(<). 


(84) 


The state-space representation of the closed-loop system is then given by: 


x(<)= [A-HC] x(t). 


(85) 
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Therefore, observing that nominal plant stability can easily be achieved, the next 
goal is to determine a systematic manner for choosing the design matrices H and G 
such that the frequency domain specifications on stability robustness and nominal 
performance (command-following, disturbance-rejection, insensitivity to sensor noise 
etc.) can be achieved. The next section shows how MBCs can be designed in a 
systematic manner so as to meet several frequency domain specifications, in addition 
to nominal closed-loop stability. 

IV.6 Linear Quadratic Gaussian/Loop Transfer Recovery (LQG/LTR) 

The LQG/LTR design method was first developed by Doyle and Stein [16]. It is 
a systematic procedure by which stability and certain frequency domain specifications 
can be met in the design of multivariable feedback control systems [15]. As discussed 
in the previous section, the LQG/LTR compensator belongs to the class of MBCs. 
The first step in the design process is design of the Target Feedback Loop (TFL). 

IV.6.1 Target Feedback Loop Design 

The TFL is designed by solving an artificial KF problem. The KF was developed 
by Kalman and Bucy, and it has found extensive use in several diverse fields [22]. 
The KF is a method by which real-time “optimal” state estimates, x(t), and outputs, 
y(<), are calculated based on past sensor measurements and applied controls actions. 
To perform this function the KF contains a model of the plant /sensor system, and 
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thus it is called a model-based observer or estimator. The mathematical form of the 
stochastic state dynamics (without control) is given by: 

*(<) = Ax(t) + L{(i); x(t) € » n ,*(t) e R m , (86) 

where ((£) is a fictitious Gaussian, zero-mean, white noise (process noise) with identity 
intensity matrix I. That is : 

E (£(<)] = 0. Vf, (87) 

cov ({(<); «r)| = E [{(!)<V)] = I«((-r). (88) 

In the measurement equation, sensor errors are modeled by additive white noise: 

y(0 = Cx(t) + 6{t) ; y(t) € e » n , (89) 

where #(t) is Gaussian, zero-mean, white noise (sensor noise) with constant intensity 
matrix pi, p > 0, independent of the process noise. That is : 

E [*(<)] = 0, Vt, (90) 

cov [<((); #(t)] = E [(((iCfr)] =^U(t-r), (91) 

and 

cov[£(l)i#(r)] = £ [{(O^r)] =0. (92) 

A KF is driven by the control vector, u(t), and the noisy measurement vector, y(£ ), 
and generates in read-time a state estimate vector, x(t) 6 J? n , and an output estimate 
vector, y(£) 6 The KF equations are : 

(93) 


x(£) = Ax(£) + Bu(f) + H [y(t) — C x(£)] 

y (<) = c*(f), 


(94) 



or 
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*(0 = [A - HC] x(<) + Hy(<) + Bu(t), 

y(0 = Cx(f), 

where H is the KF gain matrix, which is given by: 


(95) 

(96) 


H = (^C' (97) 

where E, the error covariance matrix, is determined by solving the Filter Algebraic 
Riccati Equation (FARE) given by: 


0 = AE + E A' + LL' - (i)EC'CE. (98) 

If [A,L] is stabilizable and if [A,C] is detectable, then it can be shown that there 
exists a unique nxn solution matrix E to FARE, such that E is positive definite. 

In the above equations fi and L are the design parameters of the KF problem. 
These design parameters are chosen iteratively such that the KF loop (the target loop) 
satisfies certain frequency domain specifications. Specifically, the singular values for 
the TFM, G(s)K(s), are to be large at low frequencies and low at high frequencies. 
Also, the bandwidth of the system can be set at particular values by adjusting ft. 

The process of designing the KF loop such that its maximum and minimum 
singular values (<r max and <T min ) have a predetermined “character” is termed as “loop 
shaping”. In this manner, frequency domain specifications like good command fol- 
lowing, disturbance rejection and insensitivity to noise can be met [20]. Qualitatively 
speaking the singular values of the KF loop should have the following characteristics: 
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^m t n(G KF (5)) > 1, (g g 

in the low frequency range, and, 

tr mai(S KF ( 5 )) < 1, ^ 10() 

in the low frequency range, where, 

Loop TFM: G KF (s) = C(«I- A)- J H, 

Sensitivity TFM : S KF ( s ) = [I + Gkf(s)]- 1 , and, 

Closed-Loop TFM : C KP (s) = [I + G K pfj)]" 1 G KF ( S ). 

Also, for certain performance specifications to be met the plant model may have 
to be augmented with additional dynamics (integrators etc.) at the plant input. The 
filter gam matrix H is then determined using the augmented plant model. 

IV.6.2 Actual Loop Robustness Recovery 

Having determined the KF gain matrix, H, the next step is in determining the 
control gain matrix, G, snch that the properties of the KF loop, snch as command 
following, disturbance rejection and insensitivity to sensor noise are preserved when 
using the actual loop G(s)K(s). A family of control gain matrices G„ is determined 
by solving the cheap-control LQR problem for minimum phase plants, as follows [15]: 

G|-(-)BK p , asp— »0, (101) 

where, the matrix K, is the solution to the following Control Algebraic Riccati 
Equation (CARE), and p is a design parameter 

0= -K P A - A'Kp - C'C + (-)KpBB'Kr,. 

P 


(102) 
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It can be shown that G p will always be such that : 

ReA t (A - BG P ) < 0, Vp>0 ; i = 1,2, . . . ,n, (103) 


lim y/~p G p — + W C W'W = I. 

P-> 0 

Now, the forward loop transfer function matrix of the actual loop is: 


T(s) = G(i) K p (i) = C §(5 ) B G(, [$ ! («) + B G p + H C] H, (104) 

where 

*(s) = (si -A) -1 . 

From Equation (104), the implication is that as p -*■ 0: 

C $(s)B G P [*~\s) + B G p + H C] 1 H — * C *(«) H. 

This implies that as p -» 0, the MBC loop will approximate the behavior of 
the designed target loop with the assumption that the plant is minimum phase (23). 
Therefore, if the filter loop meets the required performance specifications, then the 
MBC loop will also meet the same specifications. Hence, for minimum phase plants, 
the intuitive explanation for the operation of the LQG/LTR compensator is that 
it produces an approximate inverse of the plant G(s), and it substitutes the plant 
dynamics with the desirable dynamics of the KF loop, Gkf(«)- That is, 


Kp(s) = G 1 ( 5 ) g kf(«)- 


(105) 
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IV. 7 Guaranteed Properties of LQG/LTR Compensators 


The principal advantage of designing compensators for linear MIMO systems 
using systematic procedures lies in the fact that, for minimum phase systems, several 
good properties of the designed compensators are guaranteed. This enables the 
designer to have greater flexibility in the compensator design as several of the nominal 
properties of the LQG/LTR compensator are assured. In addition to nominal stability 
some of the guaranteed properties of the LQG/LTR compensator for minimum phase 
systems are as follows [20] : 

!•) There is a guaranteed mini mum gain margin (GAf), in each loop of the general 
MIMO system, given by 


GAT = 




(106) 


2.) There is a guaranteed minimum phase margin (PM), in each loop of the general 
MIMO system, the range of which is given by: 


-60° < PM < +60°. (107) 

The above properties offer an advantage over other design methods, such as pole 
placement or LQG, for which stability robustness is not guaranteed [21]. At this point 
it is important to emphasize that the above guaranteed properties Me true only for 
minimum phase systems, and they cannot be ensured for systems with non-minimum 
phase behavior. 
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CHAPTER V 

GAIN-SCHEDULED CONTROLLER 
DESIGN FOR VOLTAGE REGULATION 


V.l Introduction 

This chapter discusses the actual LQG/LTR compensator design for the SNPS 
using the linearized reference models. The LQG/LTR compensators thus designed 
are analyzed, via investigating the singular values of the appropriate TFMs (as 
was discussed in the previous chapter), to determine if the LQG/LTR guaranteed 
properties are achieved. Finally, the gain- scheduling of the Unear compensators to 
obtain a single nonhnear external controUer is described. 

V.2 Design and Analysis of External ControUers for the Static Space Nuclear Power 

System 

This section describes the design sequence of an LQG/LTR compensator for 
the SNPS system. For this system, the design specifications are the following: 

1. ) The designed compensator should exhibit stabiUty robustness to system param- 

eter variations and unmodeled dynamics, 

2. ) The compensator should be capable of rejecting disturbances and regulate the 

load voltage with zero steady-state error by appropriately manipulating the 

external core reactivity, 
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3. ) The compensator should be capable of slowing down the load voltage response 

of the system, thus allowing sufficient time for the movement of the control 
drums, and, 

4. ) The initial reactivity rate should not exceed 80 cents. 

The third design specification sets an upper limit on the bandwidth of the 
system. The singular values for the open-loop system models linearized at 80% and 
110% of full power are shown in Figures 44, 45, respectively. From Figures 44 and 
45, it is apparent that the cross-over frequencies (or the system bandwidth) are too 
high to allow for sufficient time for the movement of control drums. Furthermore, 
the figures indicate the existence of finite D.C. gain, causing non-zero steady-state 
error. Therefore, the linearized reference model is augmented with an integrator at 
the input. The last design specification is met by augmenting the compensator input 
by a low-pass filter which delays and reduces the magnitude of the load voltage error. 
The block diagram of the overall control system in series with the linearized reference 
model is displayed in Figure 46. In Figure 46, the reference input is set to be zero, 
and e(t) equals to 6V(t) which is the deviation of the load voltage from 100 Volts, 
e f (t) is the filtered voltage error, p ex (t) is the external reactivity is the system 

disturbance which is proportional to the change in electrical load demand, and finally, 
n(<) is the output measurement noise, which is modeled as white Gaussian noise with 
a certain variance. As it can be seen from this figure, the compensator, K( s ) is 
augmented by a filter in order to filter the load voltage error, e(t), so as to prevent 
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Figure 44. Singular Values of the Open- Loop SNPS Model Linearized at 80% of Full 
Power. 
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Figure 45. Singular Values of the Open-Loop SNPS Model Linearized at 110% of 
Full Power. 
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the initial reactivity rate from exceeding 80 cents/sec. The filter transfer function, 
Gf(s) is given by: 

Gf(s) = Cf (si - A f ) J B f (108) 

where the filter transfer function matrices, A f , B f , and C f are given in Appendix 
C. The plant, G(s), is augmented with an integrator in order to minimize the load 
voltage steady-state error. The compensator transfer function K(s) is determined 
using the step by step LQG/LTR design method described in the previous chapter. 
In Figure 46, the compensator transfer function, K(s) is given by: 

K(s) = G(sl — A + BG + HC)" l H. (109) 

The two steps of the LQG/LTR design method, which uses the augmented 
linearized reference model of the SNPS to determine K(s), are as follows: 

1.) The filter gain matrix H is fully determined by choosing a value for y and 
L, and by solving the FARE (Equation (98)). L is set to the input vector 
B aug . Initially, y is chosen to be 10® and the singular values of the target 
loop* <7‘[C(yu;I — A) H], are plotted to determine the cross-over frequency. 
If the resulting target loop does not exhibit the desired characteristics, y 
is changed and the above calculations are repeated until the singular values 
have the required specifications. In the present compensator design, y is 
chosen such that for each linear compensator designed to control the system at 
different equilibrium points, the cross-over frequency of the target loop system 
is about 0.16 rads/sec. This cross-over frequency corresponds to a rise time of 
approximately 7.0 seconds. 
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Table 5. Values of y and Corresponding Bandwidths of the Recovered Loops at 
Different Operating Points 


Model Linearized At 
50% 


fLTR (rads/sec ) 


60% 

70% 

80% 

90% 

95% 

100 % 

105% 

110 % 

120 % 

130% 

140% 

150% 


2.3700E+09 

2.6680E+09 

2.9820E+09 

3.3050E+09 

3.6000E+09 

3.7430E+09 

3.8380E-I-09 

3.9920E+09 

4.1820E+09 

4.5480E+09 

4.8990E+09 

5.2300E+09 

5.5590E+09 


1.59872E-01 

1.59892E-01 

1.59912E-01 

1.59912E-01 

1.59912E-01 

1.59912E-01 

1.59912E-01 

1.59932E-01 

1.59912E-01 

1.59952E-01 

1.59952E-01 

1.59952E-01 

1.59932E-01 


2.) Having determined the filter gain matrix H, the next step involves execution of 
the recovery procedure to satisfy the following asymptotic condition: 

lim [G aug (ja,) K tQG/LTR (ju/)] = C{jul - A) -1 H, (110) 
which implies that 

lim [ G aug (jo;)K L QG/ LT R(ju;)] = G K fO'u/). (HI) 

Hence, for p 0 the CARE, Equation (102), is solved to compute the control 
gain vector G. 

As mentioned earlier, the bandwidths of the target loops are fixed to be about 
0.16 rads/sec at all operating points. The bandwidths of the recovered loops at 
different operating points and the values of p chosen for the corresponding target 
loops are depicted in Table 5. The final gain matrices, H and G, obtained from the 
above design process are presented in Appendix D. 
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Figures 47 and 48 show the extent to which the target characteristics can 
be recovered by the LTR procedure for the 80% and 110% of full power linearized 
reference models, respectively. Since the plant transfer function matrix, G(s), is 
minimum-phase, G aug (;u>) K L q G / LXR (j'u>) of the system approximates the transfer 
function matrix of the target loop, GKp(iw), as can be observed from Figures 47 and 
48. Therefore, the actual loop has the command-following and disturbance-rejection 
properties of the target loop. It is further observed that LQG/LTR method yields 
good agreement between the singular values of the actual and the target loop for 
frequencies well beyond crossover. The singular values of GkfCM roll-off at -20 
db/dec, while those of G aug (ju;) Klqg/jjp.C/u;) eventually roll-off at -40 db/dec. 
Thus, the actual LQG/LTR loop offers some additional robustness to high frequency 
unmodeled dynamics as compared to the target loop. 

The singular values of the sensivity and the closed-loop tranfer function matrices 
can be further observed to investigate the performance and robustness properties of 
the system. These properties are achieved if: 


&maz [S(ju;)] £ 0 db. 

(112) 

&maz [C(jU>)] ^ 6.02 db) 

(113) 


where S(ju>) and C(ju> ) are the sensitivity and closed-loop transfer function matrices 
given, respectively, by: 


S(>u/) = [I + G(ju,)K(ju/)]"\ 


(114) 
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Figure 47. 


Values of the Target Loop and the Actual Recovered 
the 80% of Full Power Linearized Reference Model. 


Loop for 
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Figure 48. Singular Values of the Target Loop and the Actual Recovered Loop for 
the 110% of Full Power Linearized Reference Model. 
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and 

C(;u0 = [1+ GOu/JKOu/)]" 1 [G(jw)K(>u/)]. (115) 

The singular values of the sensitivity transfer function matrix for the 80% and 110% 
of full power linearized reference models are plotted in Figures 49 and 50. In 
Figures 49 and 50, the maximum singular values of the sensitivity transfer function 
for the 80% and 100% of full power linearized reference models are -3.333 x 10~ 7 

db and -6.457 x 10" 6 db, respectively. Hence, the condition set by Equation (112) is 
satified by the actual system. 

The singular values of the closed-loop transfer function matrix for the 80% and 
110% of full power linearized reference models are shown in Figures 51 and 52. In 
Figures 51 and 52, the maximum singular values of the closed-loop transfer function 
for the 80% and 110% of full power linearized reference models are 3.826 db and 3.434 
db, respectively. Therefore, Equation (113) is satisfied, and the designed system will 

not amplify the disturbances, or become unstable to multiplicative modeling errors 
reflected at the plant output. 

Thus, from the above analysis it is observed that the linear compensators retain 
several of the properties guaranteed by the LQG/LTR method. 

V.3 Gain-Scheduling of the External Controllers 

As it was outlined in previous chapters, linearized reference models of the SNPS 
are obtained at several operating points covering its entire operating range. A linear 
compensator is designed for each linearized reference model of the SNPS, such that 



G(jcj) K(|u)|" | (db) 







G(jcj) K(jcj)|-'[G(jcj) K(jcj) | ] (db) 
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Figure 51. 


Singular Valuts of [I + G(ju/) K(jw)]~ l (G(>u.)KO'w)l 
Full Power Linearized Reference Model. 


for the 80% of 
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when perturbed about its operating points the feedback system is nominally stable, it 
exhibits desired nominal performance, and some stability robustness. The technique 
by which the gains of the linear compensators are interpolated (scheduled) between 
these operating points in order to obtain a single nonlinear compensator that is valid 
for the entire operating range of the SNPS is known as gain-scheduling. It can be 
shown that, even though the overall gain scheduled compensator need not have any of 
the nice properties of each of the linear time-invariant compensators, the time- varying 
compensator will have guaranteed robust stability and nominal performance provided 
that [24]: 1) the gain-scheduling variable does not vary rapidly with time, and, 2) the 
small nonlinearities of the system are not captured by the scheduling variable. 

The gain vectors H and G obtained for each of the linear compensators 
(corresponding to different power levels) are valid for deviations for which the linear 
model is valid. In the gain-scheduling process these linear compensators, designed 
to control the SNPS within a finite range about its respective equilibrium points, 
are scheduled to arrive at a single compensator. This scheduled compensator can 
eliminate the transients which could result from modeling errors, while continuously 
regulating the load voltage of the SNPS. 

As the operating point about which the nonlinear reference model was linearized 
was the system electrical power, the gain-scheduling variable used in this thesis is the 
electrical power generated by the system itself. This variable is also closely related to 
the reactor thermal power, which is indicative of the system’s changing conditions. 



116 


In addition to gain vectors, the elements of the system matrices (i.e. A, B. C, L) 
are also fitted as a function of P e because the compensator uses the system dynamics 
explicitly. Therefore, the compensator state equations can be expressed as: 

z c (t) = A c (P e )z c (0 + B c (/> e ) {y(t) - y ref ) (t) ;z c (t)e» n , (11 6 ) 

Su(t) = C c (P e )z c (t), (U7) 

where z c (t) are the states of the compensator and A c ,Bc,C c are the compensator 
system matrices with elements that are functions of P e . The input to the compensator 
>s the filtered feedback loop error (the filtered load voltage deviation), and the 

compensator output is the SNPS input (the external reactivity). The compensator 
matrices are expressed as follows: 


A c (P e ) = A (P e ) - B(P e ) G(P e )~ U(P e ) C (P e ), 


( 118 ) 


Bc(-Pe) = -H(/> e ), 


( 119 ) 


C c (P e )= -G(P e ), 


( 120 ) 


where, A(P e ), B(P e ), and C(P e ) are the scheduled reference model matrices given in 

Appendix B, and G(P e ) and H (P e ) are the scheduled gain matrices given in Appendix 
E. 
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CHAPTER VI 

COMPUTER SIMULATION RESULTS 


VI. 1 Introduction 


In this chapter the main computer simulation results are presented. The mea- 
surable state and output responses of the SNPS using a linear controller are presented 
along with the control input (external reactivity) generated by the linear compensator. 
In the evaluation of the gain-scheduled compensator, the simulations include the sen- 
sor (measurement) noise effects on the load voltage of the SNPS. Finally, robustness 
of the gain-scheduled compensator is investigated through simulations which take into 
account the drifts in the system’s operating parameters, allowing investigation of the 

extent to which these variations affect the load-following capabilities of the overall 
system. 

VI-2 Evaluation of the Linear Controllers 

In this section, the performance of the linear controllers is investigated via 
simulations of the SNPS’s measurable state and the output responses to an external 
reactivity input generated by the controller to meet the demanded electrical power. As 
mentioned earlier, the measurable states of the SNPS model are the reactor neutron 
power, the core coolant inlet and exit temperatures and the total core coolant mass 
flow rate. The SNPS’s controlled output is the load voltage. 
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Table 6. Description of the Controllers Used in the Simulations 


Controller 1 

Usbs the 100% of Full Power Linear Compensator ~j 

Controller II 

Uses the 1 00% of Full Power Linear Compensator up to 
120% of Full Power, then Switches to the 120% of 
Full Power Linear Compensator. 

Controller III 

Uses the 100% of Full Power Linear Compensator untill 
80% of Full Power, then Switches to the 80% of Full 
Power Linear Compensator. 


Figures 53 through 58 show the responses of the measured states and output 
of the integrated nonlinear SNPS model to a change in the reactivity generated by 
Controller I and Controller II corresponding to a step load change from 100% to 
120% and from 120% to 150% of the full electrical power level. A description of the 
controllers used in the simulations are shown in Table 6. 

Figure 53 compares the load voltage responses of the SNPS using Controller I 
and Controller II to a step load change from 100% to 120% and from 120% to 150% of 
full power and it demonstrates that the 100% of full power linear compensator alone 

is capable of regulating the load voltage at 100 Volts with less than 0.1% relative 
error. 

Additionally, Figures 53 through 58 show that using Controller II instead of 
Controller I improves the voltage regulation and the state responses of the SNPS 
only slightly, therefore, it can be concluded that the validity range of Controller I 
(100% of full power linear compensator) extends up to and beyond 120% of the full 


power level. 
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Figure 53. Load Voltage Responses for Controller I and Controller II Corresponding 

to a Load Change from 100% to 120% and from 120% to 150% of Full 
Power. 
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Figure 54. Neutronic Power Responses for Controller I and Controller II Corre- 
sponding to a Load Change from 100% to 120% and from 120% to 150% 
of Full Power. 
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Figure 55. Core Coolant Exit Temperature Responses for Controller I and Controller 
II Corresponding to a Load Change from 100% to 120% and from 120% 
to 150% of Full Power. 
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Figure 56 . 


c ore Coolant Inlet Temperature Responses for Controller I and Con- 

Im ££? SC u " i Ch “ ge from 100% *° 120% “ d 
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Figure 57. 


Total Core Coolant Mass Flow Rate Responses for Controller I and 
Controller II Corresponding to a Load Change from 100% to 120% and 
from 120% to 150% of Full Power. 
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Figure 58. External Reactivity Generated by Controller I and Controller II Corre- 
of FVdlRjwer fr ° m 100% to 120% 111(1 from 120% to 150% 
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Figures 59 through 64 show the responses of the measured outputs of the 
integrated nonlinear SNPS model controlled by Controller I and Controller III for 
a step load change from 100% to 80% and from 80% to 50% of the full electrical 
power level. From Figure 59 it can be seen from this figure that both Controller I 
and Controller II are capable of regulating the voltage at about 100 Volts. It is 

apperant from Figures 59 through 64 that using Controller III does not improve the 
responses of SNPS over using Controller I. However, in the presence of sensor noise 
incorporated to the load voltage output of the SNPS, a gain-scheduled compensator 
can be needed for voltage regulation. 

The above simulations show that the linear compensators can successfully 
regulate the load voltage at about 100 Volts. 

VI.3 Performance Evaluation of the Gain-Scheduled Controller 

This section investigates the performance of the gain-scheduled compensator in 
the presence of sensor noise which is modeled as white Gaussian noise with a variance 
of 0.1. Figures 65 through 71 show the responses of the measured states and output 
of the nonlinear SNPS model controlled by the gain-scheduled compensator to a load 
change from 100% to 150% of the full electrical power level. The noisy measured 
load voltage output of the system in Figure 65, does not give any implications of 
the load-following of the SNPS using the gain-scheduled compensator, therefore, this 
measured load voltage is smoothened out by a filter augmented to the input of the 
compensator as mentioned in the earlier chapters, to obtain a filtered load voltage 
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Figure 59. 
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Figure 60. 


Neutroiuc Power Responses for Controller I and Controller 
sponding to a Load Change from 100% to 80% and from 80% 


III Corre- 
to 50% of 
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Figure 61. Core Coolant Exit Temperature Responses for Controller I and Controller 
III Corresponding to a Load Change Rom 100% to 80% and from 80% 
to 50% of Full Power. 
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Figure 62. Core Coolant Inlet Temperature Responses Controller I and Controller 
III Corresponding to a Load Change from 100% to 80% and from 80% 
to 50% of Full Power. 
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Figure 63. Total Core Coolant Mass Flow Rate R« nnn f n 
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Figure 64. 
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regulation error. Figure 66 shows the filtered load voltage regulation error response of 

the SNPS to a change the react, vity generated by the Gatu-Scheduled Compensator 

Corresponding to a load change from 100% to 120%. and from 120% to 150% of 

f«H power. It can be seen from th.s figure that the gain-scheduled compensator can 

regulate the load voltage with almost rero filtered load voltage regulation error desptte 
noisy sensor readings. 

Figures 72 through 77 show the responses of the measured states and output 

of the nonlinear SNPS model to an external reactivity change generated by the 

ed tiled compensator which corresponds to a load change from 100% to 80% 

end from 80% to 50% of the full electrical power level. Figure 72 shows that the 

filtered load voltage regulation error is stabilised a. about aero by the gain-schednled 

mpensator for a load change from 100% to 80% and from 80% 50% of full power. 

I. can be concluded from the simulations that the gain-scheduled compensator 

can regulate the load voltage at 100 Volts in the presence of sensor noise with aero 
steady-state regulation error. 


VI.4 Robustness of the Gain-Scheduled Controller 

Figures 78 through 89 show the responses of the measured states and output 
of the integrated nonlinear SNPS model to a change in react, vtty generated by the 
gain-scheduled compensator corresponding to a step M change from 100% to 120% 


and from 120% to 150% of the full electrical 


power level for various cases of variations 


in the fuel feedback reactivity coefficient. In these figures three 


cases are investigated: 
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Figure 66. Filtered Load Voltage Regulation Error for the Gain-Scheduled Compen- 
sator Corresponding to to a Load Change from 100% to 120% and from 
120% to 150% of Full Power. 
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Figure 67. Neutronic Power Response for the Gain-Scheduled Compensator Corre- 
sponding to a Load Change from 100% to 120% and from 120% to 150% 
of Full Power. 



Core Coolant Exit Temperature (K) 


136 



Figure 68. 


f* tke Gain-Scheduled Com- 
120% to 150% of Full | ower Cha ° 8e from W0% lo 120% and from 



Core Coolant Inlet Temperature (K) 


137 



Figure 69. Core Coolant Inlet Temperature Response for the Gain- Scheduled Com- 
pensator Corresponding to a Load Change from 100% to 120% and from 
120% to 150% of Full Power. 
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Figure 70. Total Core Coolant Mass Flow Rate Response for the Gain-Scheduled 
Compensator Corresponding to a Load Change from 100% to 120% and 
from 120% to 150% of Full Power. 
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Figure 72 . 


Filtered Load Voltage Regulation Error Response for the Gain-Scheduled 
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Figure 73. 


Neutronic Power Response for the Gain-Scheduled Compensator Corre- 
spondmg to a Load Change from 100% to 80% and from 80% to 50% of 
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Figure 74. Core Coolant Exit Temperature Response for the Gain-Seheduled Com- 

a cwee ioo% - 




Figure 75. Core Coolant Inlet Temperature Response for the Gain-Scheduled Com- 
pensator Corresponding to a Load Change from 100% to 80% and from 
80/O to 50% of Full Power. 
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Figure 76. 
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1) Constant fuel temperature coefficient (a, = -1 0 x 10~ 5 t *>\ a 

1 AU x 1U ) , 2) A negative ramp 

Wh.ch changes .he a/ from -1.0 x 10- .o -3.0 x ,0- (Negative Ramp, in I200 
seconds, 3) A pos.tive ramp which changes the a, from -1.0 x , 0 -« to zeto (Positive 

R-p, , 200 seconds. For each of these three cases, the electnca, ioad demand is 

increased from 100% to 150% of full power in two steps. In the 8rs, step, the load 

demand changes from 100% to 120% of full power, the next step i, changes from 

120% to 150% of full power. As seen from Figures 78 and Figure 79 the case which 

considers a negative ramp change in the „„ the load voltage is stabilised around 98 

Volts with a 2% error tom 100 Volts, whereas for the positive ramp case in which the 

increases linearly to rero, the load voltage becomes unstable after a certain value 
of the a f (-2.5 x 10 -6 ). 

The effects of ramp variations in the a, on the external reactivity generated by 
the gain-scheduled compensator can be seen in Figure 89. For the constant a, case, 
the external reactivity reaches a s.eady-state mine of 23 cents, whereas the external 
react. vtty generated in the negative ramp case increases in order to overcome the 
high negative fuel feedbadt reactivity resulting Rom a ramp decrease in the The 

external reactivity generated by the controller for the positive ramp case approaches 
zero, as the a f varies linearly approaching zero. 

Figures 90 through 95 show the responses of the measured states and output of 
•he nonlinear SNPS mode, to a change in reactivity generated by the gain-scheduled 
compensator to a load change from 120% to 150% of the full electrical power level 
for various cases of variations in the fuel temperature feedback reactivity coefficient. 
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Figure 78. Filtered Load Voltage Regulation Error Responses for Various Changes in 
the Fuel Temperature Feedback Reactivity Coefficient to a Load Change 
from 100% to 120% of Full Power. 
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Figure 80. Neutronic Power Responses for Various Changes in the Fuel Temperature 
Feedback Reactivity Coefficient to a Load Change from 100% to 120% 
of Full Power. 
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Figure 81. Neutronic Power Responses for Various Changes in Fuel Temperature 
Feedback Reactivity Coefficient to a Load Change from 120% to 150% 
of Full Power. 
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Figure 82. Core Coolant Exit Temperature Responses for Various Changes in the 
Fuel Temperature Feedback Reactivity Coefficient to a Load Change from 
100% to 120% of Full Power. 
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Figure 83. Core Coolant Exit Temperature Responses for Various Changes in the 

f onl TemPer t tUre FeCdback Re ^tivity Coefficient to a Load Change from 
120% to 150% of Full Power 8 
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Figure 84. Core Coolant Inlet Temperature Responses for Various Changes in the 
Fuel Temperature Feedback Reactivity Coefficient to a Load Change from 
100% to 120% of Full Power. 
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Figure 85. Core Coolant Inlet Temperature Responses for Various Changes in the 

fool ,^‘ U ? F “ dUck R«'tivity Coefficient to a Load Change from 
120% to 150% of Full Power. 
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Fire 86. Total Core Coolant Mass Flow Rate Responses for Various Changes in 
the Fuel Temperature Feedback Reactivity Coefficient to a Load Change 
from 100% to 120% of Full Power. 
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Figure 87. Total Core Coolant Mass Flow Rate Responses for Various Changes in 

r Fi ;ll emperat ^ re Feedback Reactivity Coefficient to a Load Change 
from 120% to 150% of Full Power. 8 
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Figure 89. External Reactivity Generated by the Gain-Scheduled Compensator for 
Various Changes in the Fuel Temperature Feedback Reactivity Coeffi- 
cient Corresponding to a Load Change from 120% to 150%. 
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These cases include: 1) constant ay (constant ay), and, 2) a negative step change in 
ay from its original value of —1.0 x 10 -5 to —2 x 10~ 5 (Negative Step). 

Figure 90 shows a comparison of the load voltage response of the SNPS for the 
constant cif case and negative step case. The negative step change in qj causes the 
load voltage to drop initially by about 2%, however, the system response recovers 
and finally stabilizes at about 100 Volts. The effect of a negative step change in 
ay is apparent in the neutronic power response as shown in Figure 91 where 
the neutronic power initially experiences a significant drop (-8%) from its steady- 
state value at 120% of full power level, because of high negative feedback reactivity 
resulting from changing the ay by 100%. This effect can be explained by examining 
the external reactivity generated by the compensator as shown in Figure 95. It can be 
seen from this figure that the controller responds to a negative step change in ay by 
increasing the external reactivity to overcome the high negative feedback reactivity in 
order to meet the demanded electrical power. The external reactivity generated for 
the constant ay case reaches a steady-state value of about 23 cents whereas for the 
negative step case, the external reactivity reaches a steady-state value of 45.9 cents, 
which approximately twice of that for the constant ay case. 

From the robustness evaluation analysis of the gain-scheduled compensator, it 
can be concluded that the gain-scheduled compensator can handle variations in the 
fuel temperature feedback coefficient up to a certain value. The maximum value of 
negative step change in the a y which the gain-scheduled compensator can handle in 
the presence of sensor noise has been determined to be a negative step change of 240% 
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Figure 91 . Neutromc Power Responses for Various Changes in the Fuel Temperature 
F “, ] ac n k Reactivity Coefficient to a Load Change from 120% to 150% 
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Figure 92 . 
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Figure 93. Core Coolant Inlet Temperature Responses for Various Changes in the 
Fuel Temperature Feedback Reactivity Coefficient to a Load Change from 
120% to 150% of Full Power. 
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Figure 95. External Reactivity Generated by the Gain-Scheduled Compensator for 
Various Changes in the Fuel Temperature Feedback Reactivity Coeffi- 
cient Corresponding to a Load Change from 120% to 150% of Full Power. 
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of the initial a,. The limitations are set by the fact that the maximum reactiv.ty 
insertion into the core can not exceed 80 cents. The gam-scheduled compensator can 
handle negative ramp changes in a, more successfully than positive ramp changes. 
As the a, approaches zero, as expected the system’s load voltage becomes unstable. 
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CHAPTER VII 

SUMMARY AND CONCLUSIONS 


VII. 1 Summary 

The use of shunt regulators for the load-following of the static SNPSs raises 
a number of concerns, such as the possibility of a single-point failure in the shunt 
regulators requiring reactor shutdown, or the need to deliver somewhat higher power 
than originally expected due to changing mission needs. Therefore, a back-up system 
is needed in static SNPSs to eliminate the possibility of a single-point failure in the 
shunt regulators and to improve the overall system reliability despite changes in the 
system’s operating parameters. 

Specifically, in this thesis the LQG/LTR technique, which utilizes an MBC 
structure, is used to design a compensator. LQG/LTR is a systemic MIMO control 
design methodology, which eliminates the large number of tunning parameters in ad- 
hoc methods, thus simplifying the design process. Furthermore, the LQG/LTR design 
method results in linear compensators with desirable properties, such as nominal 
stability, command following, disturbance rejection, noise insensitivity and certain 
amount of stability robustness, which though not guaranteed to carry over in the 
nonlinear domain, they appear to provide some improvements. 

The success in using the LQG/LTR method depends on the accuracy of the 
linearized models used to design it. The integrated, nonlinear SNPS model was 
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simplified to obtain a reference model which was linearized to be used in the design 
of the linear compensators. 

The LQG/LTR method is used to design linear compensators for each of the lin- 
earized reference models. The specifications for the compensator design are in terms 
of the singular values of the appropriate TFMs. Therefore, the resulting compensator 
design can be determined to be satisfactory by inspection of the appropriate singular 
values. The bandwidth of the open-loop system being the tuning variable in the com- 
pensator design is fixed at a particular value so to allow enough time for control drum 
movements. The actual performance of the linear compensators is then determined 
subjecting the integrated nonlinear SNPS model to a disturbance. This disturbance 
is a change in the demanded electric power. If the LQG/LTR compensator performs 
sufficiently well with respect to certain system specifications, such as stability of the 
system, load-following capabilities of the power system, the limitations on the initial 
reactivity insertion rate, then a satisfactory LQG/LTR compensator has been ob- 
tained. Otherwise, the design procedure is repeated so as to influence the singular 
values of the appropriate system TFMs. At this stage the bandwidth of the system 
can be altered if so desired. Also, if the situation demands, additional dynamics such 
as integrators etc. can be augmented to the compensator input and/or output. This 
is usually the case when there is a physical limitation on the system and/or there is a 
steady-state error in the output of the system. It should be pointed out that it is easy 
to determine the direction in which the design tuning parameters sure to be changed 
The lack of a large number of tuning parameters is one of the main advantages of 
the LQG/LTR design technique. It has determined that a single linear compensator 
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designed using 100% of full power linearized model was capable of regulating the load 
voltage at the desired level. 

Once the LQG/LTR compensators are designed for each of the linearized 
reference models, the gain matrices of the compensator designs are fitted to the 
scheduling variable, the electric power output, P e . These gain matrices combined with 
the scheduled model of the SNPS developed from the linearized reference models form 
the gain-scheduled compensator. The motivation for developing a gain-scheduling 
compensator is to be able to control the power system continuously within its entire 
operating regime. The performance limits of the gain-scheduled compensator were 
determined for step and ramp variations in the fuel feedback temperature coefficient. 
The cases demonstrated include 1) a negative ramp change in the fuel feedback 
reactivity coefficient, 2) a positive ramp change upto zero, and 3) a negative step 
change of 100% of the initial value of the fuel temperature feedback reactivity 
coefficient. The second case resulted in an unstable behaviour in the load voltage 
as the fuel temperature feedback coefficient approached zero. In the other two cases 
investigated, the load voltage of the SNPS has been stabilized at the desired power 
level with small relative error. Thus, the gain- scheduled compensator demonstrated 
robustness and high performance in overcoming the drifts in the system’s operating 
parameters. 

VII.2 Conclusions 

From the performance analysis of the gain- scheduled compensator the following 


conclusions cam be drawn: 
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1. ) The linearized reference models have been found to be highly accurate and 

valid for large variations from their respective equilibrium points. Therefore, 
the 100% of full power linear compensator alone could regulate the load voltage 
at the desired level for the entire operating regime of the power system. 

2. ) The gain-scheduled compensator handled large step and ramp changes in the 

fuel temperature reactivity coefficient, in the presence of sensor noise incorpo- 
rated to the load voltage of the SNPS under consideration. 

3. ) The robustness and performance analysis of the gain-scheduled compensator 

showed that the dual-loop control concept, which incorporates an adaptive 
controller with the gain-scheduled compensator as proposed earlier, can be 
replaced by the gain-scheduled compensator itself, hence eliminating the need 
for a hardware or software intensive, complex adaptive controller. 

VII.3 Implemantation Issues 

For the actual implementation of LQG/LTR compensators on small-scale digital 
systems, the finite memory, relatively slow speeds, and the expense of the hardware 
must be considered in the overall design process. Further details of these important 
considerations have been dealt with by Moroney [25]. However, some of these 
considerations can still be incorporated in the systematic LQG/LTR design procedure. 

VII.4 Recommendations 

The following are recommandations for further work on the voltage regulation 


problem in static SNPSs : 
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1. ) The integrated nonlinear SNPS model can be made more amenable to investi- 

gating variations in the system’s operating conditions, besides changes that can 
occur in the fuel temperature feedback coefficients. Such effect could include 
failures in the TE-EM pumps, and loss of radiator’s surface area. 

2. ) In this research, only sensor noise has been incorporated on the load voltage 

measurements. The performance of the gain-scheduled compensator could also 
be investigated by incorporating process and actuator noise. 

3. ) Different linear control system design techniques, such the LQ-servo, can be 

considered for the load-following of the static SNPSs. 

4. ) Since the gain-scheduled compensator is designed for an inherently stable 

system, another controller designed for an inherently unstable system can be 
designed and used for voltage regulation in cases where the fuel temperature 
feedback reactivity coefficient is positive. 
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APPENDIX A 

LINEARIZED REFERENCE SNPS MODELS 

This appendix describes the linearized reference models of the SNPS about 
twentyone operating points. 

A.l Scaling of the Linearized Reference Model States 

Scaling of the system states is done to improve the numerical characteristics 
of the A matrix, such as its condition number which is related to stiffness of the 
state equations. Three of the linearized reference model states, namely, the neutronic 
power, and the two delay-neutron precursor concentrations are scaled by the steady- 
state value of the neutronic power at each equilibrium point. 

A.2 The Linearized Reference Model Descriptions 

The linearized reference models are described by the operating power level 
(steady state) about which the nonlinear reference model was linearized. For example: 
150% model means that the model was obtained by linearizing the nonlinear reference 
model while operating at 150% of full power. The B, and L vectors are constant for 
^ the linearized reference models, given by: 

B = [2.5e6 0000000 0], (121) 

and, 

L = [0 000000 — 4.067e— 6 0], (122) 

and D is the zero vector. The A, C and F matrices for all the linearized models are 
given in Tables 7 through 27. 



liildr 7. I Ik* System Matrices (A, O', and F) of I hr 50% of Full Power Linearized Model 
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Table 8. The System Matrices (A, C, and F) of the 55% of FHill Power Linearized Model 
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Table* 9. I lie* System Matrices (A, C, and F) of the 60% of Full Power Linearized Model 
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Table 10. The System Matrices (A, C, and P) of the 65% of Pull Power Linearized Model 
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Table 13. The System Matrices (A, C, and F) of the 80% of F\ill Power Linearized Model 
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Table M. The System Matrices (A, C, and F) of the 85% of Full Power Linearized Model 
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Table 15. The System Matrices (A, C, and P) of the 90% of Pull Power Linearized Model 
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Table 16. The System Matrices (A, C, and F) of the 95% of Full Power Linearized Model 
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Table 17. The System Matrices (A, C, and F) of the 100% of Full Power Linearized Model 
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Table 19. The System Matrices (A, C, and F) of Hie 110% of Full Power Linearized Model 
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Table 20. The System Matrices (A, C, and F) of the 115% of Full Power Linearized Model 
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Table 21. The System Matrices (A, C, and F) of the 120% of Full Power Linearized Model 
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Tnble 22. The System Matrices (A, C, and P) of the 125% of Full Power Linearized Model 
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Jahle 23. I lie System Matrices (A, C, and F) of the 130% of Full Power Linearized Model 
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Table 2-1. The System Matrices (A, C, and P) of the 135% of Full Power Linearized Model 
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Table 25 The System Matrices (A, C, ami F) of the 140% of Full Power Linearized Model 
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Table 2(5. The System Matrices (A, C, and F) of the 145% of Full Power Linearized Model 


195 


tu 

8 

© 

o 

8 

+ 

LU 


8 8 

LU LQ 

8 8 

o 5 

© © 

8 § 

LU LU 


8 8 

LU LU 

§ 8 

© o 

0 © 

1 ? 

tu LU 


CM 

O o 


8 8 s 


U 3 tU LU U 3 

88 ^ 

© © o — 

0 0 —'™ 


9 f 9 SI 

LU pi LU nj 


+ 

LU 


© © 

1 1 

LU LU 


r- 

os 


Os 

o 

<N 


Os 


o © 


8 
+ 
LU 

8 8 8 1 

© © © — ' 

III? 

I 9 ! 


fN 

O 

LU 

3 

vo 


- © ™ - 

f?n 

LU LU tU LU 

o ^ 


LU LU 

© i 


LU LU 


© © 


9 ? o 

s I « 

<N S g 


w ^ 
*U in 


? 5 5 § 9 9 5 

§ QJ U U U U U m 

P O oo ® n o 

g g oo s: rn S 

© © © 


OO , 

« r n 
© °^ <N © 
- °? Vd g 


fff?5§888 


3 

psi 


© © 


oo 

— • m 

rs os’ 


© © © © 


??5888888 

2 ||| lulu luluS 

^ © ^ © © © © o ^ 


<si £} 

© © 


i £ 1 1 i 


00 

id m 

<n <n 


© o o 


§ 

© 


w 

© 


w 

© 


LU I 


Iff is 
g I 3 1 1 


^ ^ <N 

© «n m 
V — ‘ V£> 


LU 


© © 


CU w LU 
© 


L©| ^ n n ^ tn ^ 


oo 


© 

Os 


© 

i 

CU 

fD 

<N 

© 

VC 


II 

U. 


- 

oo'db^od^b 

u 

— fSmrriDvcr-ooos 





Table 27. The System Matrices (A, C, and F) of the 150% of Full Power Linearized Model 
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APPENDIX B 

LISTING OF THE FITS FOR THE 
SCHEDULED REFERENCE SNPS MODEL 


B 1 Fits Used for the Scheduled Reference Model 

CMDc T it f< ? 0wins pa « es Iist the fits in the scheduled reference model of the 
SNPS. The fits are used to model the elements of the system matrices A, B, and C. 
The scheduling variable used is the electric power produced by the SNPS, P e 


c 

c THIS IS THE SEGMENT OF A FORTRAN PROGRAM THAT FITS THE ELEMENTS 
c OF THE LINEARIZED REFERENCE MODELS OF SNPS AS A FUNCTION OF 
c EECTRIC POWER OUTPUT OF THE SYSTEM, 
c 
c 

c THE VARIABLES USED IN THIS SEGMENT ARE: 

c am (9. 9) 'A' MATRIX ELEMENTS - STATE MATRIX 

c bm(9) 'B' VECTOR, INPUT VECTOR 

c c(9) ’C’ VECTOR, OUTPUT VECTOR 

c x SCHEDULING VARIABLE SCALED BY THE 

c NOMINAL ELECTRIC POWER PRODUCED 

c rparrm(21) VECTOR, VECTOR OF PARAMETERS 

c EXPRESSED IN TERMS OF THE STATES 

c 

c the beta(i) - 2 group 
20 amkc=l . OdO 

rpamn(2) *0 . 001614d0 
rparnn(3)=0 . 004809d0 
c the neutron generation time 
rparrm(4) *4 . 0d-07 
c the lamibda(i) - 2 group 
rparrm(5)=0 . 025853d0 
rparrm(6) *0 . 246 ldO 
c the total beta 

rparrm(7)=rparrm(2)+rparrm(3) 
c the T-H eqn constants 
rpamn(8) =amf *cpf 
if (x . le . 1 . OdO)then 

rparrm(9) *2 . 74068196d0-2 . 49397632d0*x 



*+0 . 568972018d0*x**2 
else 

rparnn(9) =2 . 790241234d0-2 . 7782255246d0*x 
*+0 . 8072848472d0*x**2 
end if 

rparrra( 10) =6 . 581062334d0-6 . 668101788d0*x 
*+1.9627778d0*x**2 
if (x .le.l . 0d0)then 
rparnn( 11 ) =-0 . 00996d0+0 . 37388d0*x 
else 

rparrm(ll)=0 . 1489629021d0+0 . 294024646d0*x 
*-0 . 0775379499d0*x**2 
end if 

if (x . le . 1 . OdO) then 
rparrm(12)=6 . 8277d0-0 .41033d0*x 
else 

rparrm(12) =6 . 69599893d0-0 . 225052886d0*x 
*-0.0519716859d0*x**2 
end if 

if (x . le . 1 . OdO) then 
rparrm(13)=9 . 6181d0-0 . 57802d0*x 

else 

rparrm( 13) «9 . 43263192d0-0 . 3170429624d0*x 
*-0 .073207497894d0*x**2 
end if 

if (x . le . 1 . OdO) then 

rparrm(l4)=10 .63032195d0-l .34773111d0*x 
*+0.3128847781d0*x**2 
else 

rparrm(14) =10 . 242d0-0 . 64923d0*x 
end if 

if (x . le . 1 . OdO) then 

rparnn( 15) =76 . 5558133d0-9 . 7058838674d0*x 
*+2 .25329290d0« , x**2 
else 

rparrm( 15 ) =73 . 761d0-4 . 6755d0*x 
end if 

if (x . le . 1 . OdO) then 

rparrm( 16) =298 . 8763343d0- 192 . 28559175d0*x 
*+68. 849275 13d0*x**2 
else 

rparrm(16) =251 . 36611922d0-92 . 256681688d0*x 
*+15 . 95455235d0*x**2 
end if 

if (x. le.l .0d0)then 

rparxm( 17) =-0 . 016985d0+0 . 63759d0*x 



rparrm(17) =0 . 2540285113d0+0 . 5014043415d0*x- 
*0 . 132226566d0*x**2 
end if 

rparrm( 1 8 ) =v*cur i 
rparnn(l8)=-4 . 0713d5+602 . 28dO*x 
rparrm(18)=0 . 56940d6*x 
rparrm(l9)=l . OdO 
rparrm ( 20 ) =t s ink 
rparrm(22)=amte*cpte 
additional rparrm for Tauin 
rparrm (23) =1 . OdO 
go to 24 
24 continue 

am(l ,1)=-1 . 0d0*(rparrm(2)+rparrm(3))/rparrm(4) 

am(l ,2)=rparrm(5) 

am(l ,3) *rparrm(6) 

am(l ,4) =pnl*alft/rparrm(4) 

am(l ,5)=0 .OdO 

am(l ,6)=0 .OdO 

am(l ,7)=0.0d0 

am( 1 ,8) =0 . OdO 

am(l,9)=0.0d0 

am(2 , 1) =(rparrm(2) ) /rparrm(4) 

am (2 , 2 ) * ( - 1 . OdO ) *rparrm ( 5 ) 

am(2,3)=0.0d0 

am (2 ,4)=0 .OdO 

am(2,5)=0.0d0 

am(2,6)=0.0d0 

am(2,7)=0.0d0 

am(2,8)=0.0d0 

am(2,9)*0.0d0 

am(3,l)=(rparrm(3) )/rparrm(4) 

am(3,2)=0.0d0 

am(3 ,3)=(-l -0d0)*rparrm(6) 

am(3,4)=0.0d0 

am(3,5)=0.0d0 

am(3,6)*0 .OdO 

am(3,7)=0.0d0 

am (3 ,8) =0 .0d0 

am(3 ,9)=0 .0d0 

am(4,l)=l .0d0/rparrm(8) 

am(4,2)=0.0d0 

am(4,3)=0.0d0 

am(4,4)=-l .0d0/rparrm(9) 

am(4 , 5 ) =0 . 5d0/rpaxrm(9) 
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am(4.6)=0.0d0 

am(4 ,7) =0 . 5d0/rparrm(9) 

am(4,8)=0.0d0 

am(4,9)=0.0d0 

am (5 , 1) =0 . OdO 

am(5 ,2)=0 .OdO 

am(5,3)=0.0d0 

am(5 ,4) =1 . 0d0/rparrm(10) 

am( 5 , 5 ) = ( - 1 . OdO ) * ( 1 . OdO*rparrm( 1 1 ) +0 . 5dO/rparrm( 10 ) ) 
am(5,6)=0.0d0 

am(5 ,7)=-0.5d0/rparrm(10)+rparnn(ll) 

am(5,8)=0.0d0 

am(5,9)=0.0d0 

am(6 , 1 ) =0 . OdO 

am(6,2)=0.0d0 

am(6,3)=0.0d0 

am(6,4)=0.0d0 

am(6 , 5 ) =rparrm( 17) -0 . 5d0/rparnn( 12) 

am (6 ,6) *-l . OdO* (rparrm(17)+0 . 5d0/rparrm(12) ) 

am(6,7)=0.0d0 

am(6,8)=l .0d0/rparrm(l2) 

am(6,9)=0.0d0 

am(7,l)=0.0d0 

am(7,2)*0.0d0 

am(7,3)=0.0d0 

am(7,4)=0.0d0 

am(7,5)=0.0d0 

am(7,6)=l .0d0/rpamn(23) 

am(7,7)«-l . 0d0/rparrm(23) 

am(7,8)=0.0d0 

am(7,9)=0.0d0 

am(8,l)=0.0d0 

am(8,2)=0.0d0 

am(8,3)=0.0d0 

am(8,4)=0.0d0 

am ( 8 , 5 ) =0 . 5d0 /rparrm ( 13 ) 

am(8 ,6)=0 . 5d0/rparrm(13) 

am(8,7)*0.0d0 

am(8,8)=-l .0d0*(l .0d0/rparrm(13)+l .OdO/rparrm(14)) 

am(8,9)*l .OdO /rparrm (14) 

am(9,l)=0.0d0 

am(9,2)=0.0d0 

am(9,3)=0.0d0 

am(9,4)=0.0d0 

am(9,5)*0.0d0 

am(9,6)=0.0d0 



Z.U1 


am(9.7)=0.0d0 

am(9 ,8) =1 . OdO/rparrnt(15) 

am(9 , 9) =- 1 . OdO* (1 . OdO/rparrm( 15) + 1 . 0d0/rparnn(16) ) 

bm ( 1 ) =pn/ rparrm ( 4 ) 

bm(2)=0 .OdO 

bm(3)=0.0d0 

bm(4)=0 . OdO 

bm(5)=0 .OdO 

bm(6)=0 .OdO 

bm(7)=0.0d0 

bm(8)=0.0d0 

bm(9)=0 .OdO 

rparrm ( 1 ) =rhot em 

rparrm ( 2 1 ) =rhoexm 

c(l)=O.OdO 
c(2)=0.0d0 
c(3)=0.0d0 
c (4) =0 . OdO 
if (x . le . 0 . 95d0) then 

c(5)=0 . 24595 128d0+0 . 129600564d0*x-0 . 04953809777d0*x**2 

else 
end if 

if ((x.gt .0.95d0) .and. (x.le.l.05d0)) then 

c(5)=0.3246625d0 

else 

end if 

if (x . gt . 1 . 05d0) then 

c (5) =0 . 28725935d0+0 . 046055781d0*x-0 . 0095683919d0*x**2 

else 

end if 

c(6)=0.0d0 

c(7)=0.0d0 

if (x.le.l.OdO)then 

c (8)=-0 . 164194325d0-0.0266547088d0*x-0.0153233547d0*x**2 

else 

c (8) =-0 . 201684183d0+0 . 0440687 154d0*x-0 . 0487757874d0*x**2 

end if 

c(9)=0.0d0 
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APPENDIX C 

DESCRIPTION OF THE NOISE REJECTION FILTER 


C.l The description of the system matrices of the noise rejection filter, Gf(s) 

The filter matrices, Af, Bf, and Cf are scalars and have the following values: 


Af = [—0.50], 

(123) 

B f = [0.25], 

(124) 


and. 


C f = [0.20]. 


(125) 
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APPENDIX D 

THE GAIN VECTORS OF THE LINEAR COMPENSATORS 

D.l The Gain Matrices of the Linear Compensators 

In this appendix the numerical values of the gain vectors H and G for thirteen 
different linear compensators are given in Tables 28 and 29, respectively. 



Table 28. Filter Gain Vector, H, Obtained for the Different Linear Compensators at Different Power Levels 
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APPENDIX E 

LISTING OF THE FITS FOR THE 
GAIN-SCHEDULED COMPENSATOR 


E l Fits Used for the Scheduled Reference Model 

The following pages list the fits used to implement gain-scheduled compensator. 
I he fits are used to model the elements of the system matrices A, B, and C of the 
reference model. The gain matrices H and G are also fitted. The scheduling variable 
used is the electric power produced by the SNPS, P e . 


c 

C THIS IS THE SEGMENT OF A FORTRAN PROGRAM THAT FITS THE 
c ELEMENTS OF THE LINEAR COMPENSATORS AS A FUNCTION OF 
c THE ELECTRIC POWER PRODUCED BY THE SYSTEM, 
c 

c THE VARIABLES USED IN THIS SEGMENT ARE: 
c r(132) — ELEMENTS OF THE LINEAR COMPENSATOR 

c h(10) FILTER GAIN VECTOR EXPRESSED IN 

c TERMS OF THE r(.)'s 

c g(lO) CONTROL GAIN VECTOR EXPRESSED IN 

c TERMS OF THE r(.)'s 

c x SCHEDULING VARIABLE SCALED WITH THE 

c NOMINAL ELECTRIC POWER OUTPUT OF 

c THE SYSTEM 


if (x.le.l. OdO) then 

r (1) =-241 . 77656423d0-125 . 82813557d0*x- 
*113. 21292827dO*x**2 
els* 

r ( 1 ) =- 1 33 . 76d0-347 . 43d0*x 
end if 

if (x.le.l .OdO) then 

r ( 2) =-0 . 0174896478d0+8 . 75806245d-3*x-0 . 037150464d0*x**2 

else 

r (2) =8 . 70330778d-4-0 . 025337685d0*x-0 . 021621332d0*x**2 
end if 

r(3)=0.0d0 
r(4)=0.0d0 
if (x.le.l. OdO)then 



r (5) =-8 . 72496835d0+l . 70332824d0*x-7 . 3681 171 2dO*x**2 
else 

r (5) =0 . 50195175d0-16 . 363583d0*x+l . 53105584d0*x**2 
end if 

if (x . le . 1 . OdO ) then 

r (6) =-2409 . 4821706d0- 1419 . 0829657d0*x+595 . 14332435d0*x**2 
else 

r(6)=-2996.263463d0-238 .93897943d0*x+10.475350112d0*x**2 

end if 

if (x.le.l.OdO)then 

r (7) =1 . 5441046d0-0 . 90103442d0*x+0 . 28934191d0*x**2 
else 

r (7) =1 • 58272175d0-0 . 98353868d0*x+0 . 335526696d0*x**2 
end if 

if (x.le.l. 0d0)then 

r (8) =3 . 64205606d0-14 . 711441275d0*x+8 . 434039 14d0*x** 2 
else 

r (8) =- 1 1 . 469d0+8 . 8356d0*x 
end if 

if (x.le.l.OdO)then 

r(9) =1641 . 9713094d0+254 . 15676338d0*x+161 . 60976489d0*x**2 
else 

r (9) =2017 . 0039659d0-450 .37l32179d0*x+493 . 12719647d0*x**2 
end if 

r(l0)=2 . 8609797858d0-l . 580192198d0*x+0 . 5035063073d0*x=*2 

r(ll)=0.0d0 

if (x.le.l. OdO) then 

r(l2)=2 .90246324d-5-2.096244951d-5*x+8.0481963644d-6*x**2 

else 

r ( 12) =2 .7155547 d-5-1 . 46730506d-5*x+3 . 67507453d-6*x**2 
end if 

r(l3)=2.5d6 
r(l4)=-1.60575d4 
r(l5)=0.025853d0 
r(16)=0.2461d0 
r(l7)=-25 .OdO 
if (x.le.l .OdO)then 

r ( 18) =-0 . 0027886485d0+0 . 0017706142dO*x-6 . 6748138d-4*x**2 
else 

r ( 18 ) =-0 . 0027480766d0+0 . 00136594377d0*x-3 . 02178662d-4*x**2 
end if 

r(19)=0 .OdO 
r(20)=0.0d0 
if (x.le.l .0d0)then 

r(21)=0 .0017808008d0-0 . 0013775201d0*x+6.6562939d-4*x**2 

else 



r(21)=0.0016410608dO-8 . 67344804d-4*x+2 . 9781543d-4*x**2 

end if 

r(22)=0.0d0 

r(23)=0.0d0 

if (x.le.l .OdO)then 

r(24)=0.010630065d0-0.0093263183d0*x+0.0038834858d0*x**2 

else 

r(24) =0 . 00895748175d0-0 . 0048373576d0*x+0 . 0010731668426*x**2 

end if 

r(25)=0.0d0 

r(26)=4035.0d0 

r(27)=-0. 025853d0 

r(28)=0.0d0 

r(29)=0.0d0 

if (x . le . 1 . OdO) then 

r(30)=-114 . 86772685d0+75 . 314375616d0*x-27 . 599383096d0*x**2 
else 

r(30)=-113 . 64225635d0+59 . 835062282d0*x-13 . 27536245dO*x**2 

end if 

r(31)=0.0d0 

r(32)=0.0d0 

if (x . le . 1 . OdO) then 

r (33) =73 . 20809844d0-57 . 7014691d0*x+27 . 081 150244d0*x**2 

else 

r(33)=67 . 509081587d0-36 .989723269d0*x+12 . 156680824d0*x**2 

end if 

r(34)=0.0d0 

r(35)=0.0d0 

if (x . le . 1 . OdO) then 

r (36) =437 . 1947937d0-390 .79382385d0*x+160 . 27328293d0*x**2 
else 

r (36) =370 .0056601d0-210 . 10759902d0*x+46.938195816d0*x**2 

end if 

r(37)=0.0d0 

r(38)=1.20225d4 

r(39)=0.0d0 

r(40)=-0.2461d0 

r(41)=0.0d0 

if (x . le . 1 . OdO) then 

r(42)=-117.57496196d0+74.649282d0*x-27.7788018d0*x**2 

else 

r(42)=-116.40011835d0+58.8671043d0*x-13.118211847*x**2 

end if 

r(43)=0.0d0 

r(44)=0.0d0 

if (x . le . 1 . OdO) then 
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r (45) =75. 056071961 - 57 . 99157055 d 0 *x+ 27 . 774862583 d 0 *x **2 

else 

r (45 >=69.461 506179d0-37 . 196699811d0*x+12 .68228478d0*x**2 
end if 

r(46)=0.0d0 
r(47)=0.0d0 
if (x.le.l.OdO)then 

r (48) =448 . 17466022d0-393 . 28013148d0*x+162 . 70187l23d0*x**2 
else 

r (48) =379 . 2609485d0-207 . 8651728d0*x+46 . 444946921d0*x**2 
end if 

r(49)=0.0d0 

r(50)=21 . 7l4d0+30 . 427d0*x 

r(51)=0.0d0 

r(52)=O.OdO 

r(53) 0.446 26946 2 76d0+0. 16150150049d0*x-0.93686258384d0*x**2 

if (x.le.l . OdO)then 

r(54)=0 .068628583934d0-0. 29559371258d0*x+0.62931944246d0*x**2 

else 

r(54)=-0 . 72194d0+l . 1175d0*x 
end if 
r(55)=0.0d0 
if (x.le.l. OdO ) then 

r (56) =0 . 30417768093d0-0 . 29300248518d0*x+0 . 59764210265d0*x**2 
else 

r ( 56) =-0. 17279099907 d0+0.58101010423d0*x+0. 19858729212d0*x**2 

end if 

if (x . le . 0 . 95d0) then 

r (57) =0 . 14873473368d0-0 . 02022433039d0*x+0 . 0015317469082d0*x**2 

else 
end if 

if ((x.gt.0.95d0) .and. (x.le . 1 .05d0) )then 
r (57)=0 . 13116d0 
else 
end if 

if (x . gt . 1 . 05d0) then 

r(57)=0.21599116507d0-0.14582842787d0*x+0.061950756343d0*x**2 

else 

end if 

r(58)=0.0d0 

r(59)=0.0d0 

if (x.le.l .OdO) then 

r(60)=0.93871276585d0-0.38153433085d0*x+0.077841340437d0*x**2 

else 

r( 60 )=l .0782397381d0-0.62377925728d0*x+0. 18395982702d0*x**2 

end if 
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r(61)=0.0d0 
r(62)=0.0d0 
r(63)=0 .OdO 
r(64)=0.0d0 
if (x . le . 1 . OdO) then 

r(65)=0 . 267676544d0-0 . 257528202d0*x+0 . 52630175255d0*x**2 
else 

r(65)=-0 . 1523653644d0+0 . 5121013797d0*x+0 . 1749441954d0*x**2 
end if 

if (x . le . 1 . OdO) then 

r(66)=-0 . 316139159d0-0 . 17384025d0*x-0 .316644742d0«‘x**2 
else 

r(66)=-0 . 22123d0-0 . 58386d0*x 
end if 
r(67)=0 .OdO 
if (x.le.l. OdO) then 

r (68) =-0 . 11626056742d0+0 . 426381303d0*x-0 . 2126024288d0*x**2 

else 

r(68)=0 . 22488813863d0+0 . 038378732074d0*x-0 . 1651654236d0*x**2 
end if 

if (x . le . 1 . OdO) then 

r(69)»0. 109587181d0-0.029687181d0*x+0 . 0295689 1044d0*x**2 

else 

r (69)*0 . U631705d0-0 . 0371580942d0*x+0 .030807431323d0*x**2 

end if 

r(70)*0.0d0 

r(71)=0.0d0 

if (x . le . 0 . 95dO) then 

r (72 ) =0 . 64918377d0-0 . 22570807d0*x+0 . 1059844276d0*x**2 

else 

end if 

if ((x.gt .0.95d0) .and. (x.le.l .05d0)) then 

r(72)=0.531d0 

else 

end if 

if (x . gt . 1 . 05 dO) then 

r(72)=0 ,60301764d0-0 . 108657998d0*x+0 . 0397866829d0*x**2 

else 

end if 

r(73)»0.0d0 

r(74)=0.0d0 

r(75)=0 .OdO 

r(76)=0.0d0 

r(77)*0.0d0 

if (x.le.l .OdO) then 

r (78) *-0 . 077828932d0+0 . 352743204d0*x+0 . 140932428d0*x**2 
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else 

r (78) =0 . 06907 6382d0+0 . 4755497775d0*x-0 . 129550456d0*x**2 
end if 

if (x.le.l .0d0)then 

r(79)=-0. 102865 522d0-0. 51287762d0*x-0 .085853237d0*x**2 
else 

r(79)=-0 . 3284062266d0-0 .5042657899d0*x+0 . 1315858366d0*x**2 
end if 
r(80)=0.0d0 
if (x.le.l. 0d0)then 

r (81 ) =0 . 17456649d0+0 . 0756895688d0*x-0 . 01 2086411d0*x**2 
else 

r(8l)=0.2252048896d0-0.00738815938d0*x+0.02084031699d0*x**2 

end if 

r(82)=0.0d0 

r(83)=0.0d0 

if (x . le . 0 . 95 dO) then 

r(84)=0 . 167800199d0+0 .38669139d0*x-0 . 15720338d0*x**2 

else 

end if 

if ((x.gt .0.95d0) .and. (x.le.l .05d0)) then 

r(84)=-0 . 61282476857d0+l . 93163028d0*x-0 . 917917005d0*x**2 

else 

end if 

if (x .gt . 1 . 05d0)then 

r(84) *0 . 3896305d0+0 . 0105961918d0*x+0 . 0023816599d0*x**2 

else 

end if 

r(85)=0.0d0 

r(86)=0.Cd0 

r(87)=0.0d0 

r(88)=0.0d0 

r(89)=0.0d0 

if (x.le.l .0d0)then 

r(90)=-0.030360992d0-0.13910675d0*x+0.050487252d0*x**2 

else 

r(90)=-0.0996051345d0-0 .0236785024d0*x+0.00407918287d0*x**2 

end if 

r(91)=l .OdO 

r(92)=-1.0d0 

if (x . le . 1 . 0d0)then 

r ( 93 ) =0 . 0253887348d0+0 . 060873223d0*x-0 . 0108383838d0*x**2 
else 

r (93) *0 . 068620443d0-0 . 00979847597d0*x+0 . 0170355203d0*x**2 

end if 

r(94)=0.0d0 
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r(95)=0.0d0 

if (x . le . 0 . 95d0) then 

r(96)=0 . 147999563d0+0 . 36246270d0*x-0 . 14675276d0*x**2 

else 

end if 

if ( (x.gt .0 .95dO) .and. (x . le . 1 .05d0) )then 

r (96) =-0 . 579140314d0+l . 803128994d0*x-0 . 8570189753d0*x**2 

else 

end if 

if (x.gt . 1 . 05d0)then 

r (96)=0 . 354150196d0+0 .01393371122d0*x+4. 549830352d-4*x**2 

else 

end if 

r(97)=0.0d0 

r(98)=0.0d0 

r(99)=0.0d0 

r(l00)=0.0d0 

r(l01)=0.0d0 

if (x . le . 1 . 0d0)then 

r(l02)=0.018010301d0-0.0250815518d0*x+0.0089096534d0*x**2 

else 

r ( 102) *0 . 00426963 189d0-0 . 0029206475d0*x+4 . 34904379d-4*x**2 
end if 

if (x . le . 1 . OdO) then 
r ( 103) *0 . 051876d0+0 . 0034231d0*x 

else 

r ( 103) =0 . 0531 105105d0+0 . 00154143379d0*x+6 . 42869498d-4*x**2 
end if 

r(l04)=0.0d0 

if (x . le . 0 . 95 dO) then 

r ( 105) ®-0 . 174619597d0-0 . 0107835062d0*x+0 . 00435952586d0*x**2 

else 

end if 

if ((x.gt .0.95d0) .and. (x.le.l.05d0))then 

r(l05)=-0.180875d0 

else 

end if 

if (x . gt . 1 . 05d0) then 

r ( 105) *-0 . 1676871 1 1 ldO-O . 02170685834d0*x+0 . 00874028692d0*x**2 

else 

end if 

if (x.le.l .OdO)then 

r ( 106) =0 . 093853593d0+0 . 0128137098d0*x-0 . 0024473029d0*x**2 
else 

r( 106) =0 . 096919d0+0 . 0073056d0*x 
end if 
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r(l07)=0 . OdO 

if (x . le . 0 . 95d0) then 

r(l08)=0 . 138998897 d0+0 . 034342255d0*x-0 . 00959476153d0*x**2 

else 

end if 

if ( (x .gt .0 .95d0) . and . (x . le . 1 . 05d0) )then 

r (108) =-0 . 0466983878d0+0 . 399036999d0*x-0 . 187599509d0*x**2 

else 

end if 

if (x . gt . 1 . 05d0) then 

r(l08)=0 . 171449193d0-0 .0124447735d0*x+0 .006418066221d0*x**2 

else 

end if 

r(l09)=0.0d0 
r (110) =0. OdO 
r(lll)=0.0d0 
r(H2)=0.0d0 
r(ll3)=0.0d0 
if (x.le.l .OdO) then 

r ( 1 14) =-0 . 00167 1924d0-0 . 00258061887d0*x+7 . 01570863d-4*x**2 
else 

r(ll4)=-0.0026799402d0-9.50021147d-4*x+7.90890137d-5*x**2 
end if 

r(115)=0.0d0 

r(116)=0.0d0 

if (x.le.l .OdO) then 

r (117) =0 . 014349d0+0 . 0023788d0*x 

else 

r(ll7) =0 . 0154564814d0+6 . 29356966d-4*x+6 .47303304d-4*x»*2 
end if 

if (x.le.l .OdO) then 

r(118)=-0.016038793312d0-0.005217989d0*x+0.00107811734d0*x**2 

else 

r(ll8)=-0.017152d0-0.0030301d0*x 
end if 

r(ll9)=0.0d0 
if (x.le.l .OdO) then 

r ( 120) =0 . 00738678741d0+0 . 0043190172944d0*x~7 . 88904557d-4*x**2 
else 

r(120)=0.00897486168d0+0.002263234d0*x-2.88498892d-4*x**2 
end if 

if (x . le . 1 . OdO) then 

r(121)=-241 .77656423d0-125 .82813557d0*x- 
*113.21292827d0*x**2 
else 

r(l21)=-133.76d0-347.43d0*x 
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end if 

if (x.le.l. OdO) then 

r(122)=-0 . 0174896478d0+8 . 75806245d-3*x-0 .037150464d0*x**2 
else 

r(122) =8 . 70330778d-4-0 . 025337685d0*x-0 . 021621332d0*x**2 
end if 

r(l23)=0.0d0 
r (124) =0. OdO 
if (x.le.l .OdO)then 

r( 125) =-8.72496835d0+1.70332824d0*x-7. 3681 1712dO*x**2 
else 

r(125)=0 . 50195175d0-16 . 363583d0*x+l . 53105584d0*x**2 
end if 

if (x . le . 1 . OdO) then 

r (126) *-2409 . 4821706d0-1419 . 0829657d0«-x+595 . 14332435d0*x**2 

else 

r(l26) *-2996 . 263463d0-238 . 93897943d0*x+10 . 475350112d0*x**2 
end if 

if (x.le.l .OdO) then 

r(l27) *1 . 5441046d0-0 . 90103442d0*x+0 . 2893419 ld0*x««* 2 

else 

r(l27)*l . 58272175d0-0 . 98353868d0*x+0 . 335526696d0*x**2 
end if 

if (x . le . 1 . OdO) then 

r (128) *3 . 64205606d0- 14.71 1441275d0*x+8 . 434039 14d0*x** 2 
else 

r ( 128) *-l 1 . 469d0+8 . 8356d0*x 
end if 

if (x.le.l .OdO) then 

r(129)=1641 .9713094d0+254. 15676338d0*x+161 .60976489d0*x**2 
else 

r (129) =2017 . 0039659d0-450 . 37132179d0*x+493 . 12719647d0*x**2 
end if 

r (130) =2 . 8609797858d0-l . 580192198d0*x+0 . 5035063073d0*x**2 
r(131)*0.0d0 
r(132)*0.0d0 
c 

c FILTER GAIN VECTOR 
c 

h(l)*r(12) 
h(2)*r (24) 
h(3)*r (36) 
h(4)=r(48) 
h(5)*r(60) 
h(6)=r(72) 
h(7)=r(84) 



h(8)=r(96) 

h(9)=r(l08) 

h(l0)=r(l20) 

CONTROL GAIN VECTOR 

g(l)=r(l) 

g(2)=r(2) 

g(3)=r(3) 

g(4)=r(4) 

g(5)=r(5) 

g(6)=r(6) 

g(7)=r(7) 

g(8)=r(8) 

g(9)*r(9) 

g(10)*r(10) 
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